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Abstract 

We consider the dynamics in a neighborhood of an elliptic equilibrium point with a 
Diophantine frequency of a symplectic real analytic vector field and we prove the following 
result of effective stability. Generically, both in a topological and measure-theoretical 
sense, any solution starting sufficiently close to the equilibrium point remains close to it 
for an interval of time which is doubly exponentially large with respect to the inverse 
of the distance to the equilibrium point. We actually prove a more general statement 
assuming the frequency is only non-resonant. This improves previous results where much 
stronger non-generic assumptions were required. 
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1 Introduction 

The aim of this paper is to study the effective stability of elliptic equilibrium points in 
Hamiltonian systems. Our main result will be that the flow of a real analytic Hamiltonian 
H in n degrees of freedom having a Diophantine equilibrium point at the origin is doubly 
exponentially stable at the origin under an open and dense condition of full Lebesgue measure 
which only involves the part of the power expansion of H that contains the terms of degree 
between 3 and This result will be derived from a more general effective stability result 

for non-resonant elliptic equilibrium points. Before stating the exact results, let us start by 
describing the general setting. 

1.1 Stability of elliptic equilibrium points 

We consider a symplectic manifold (M, H) of dimension 2n, n € N, where H is an everywhere 
non-degenerate closed 2-form, a smooth symplectic vector field X on M (meaning that the 
one-form ix^ is closed, or, equivalently, that the Lie derivative Cx^ vanishes identically) 
and an equilibrium point p* £ M, that is X{p*) = 0. We are interested in studying whether 
p* is stable in the following sense (in the sense of Lyapounov): given any neighborhood U of 
p*, there exists a smaller neighborhood V of p* such that for any point po G V, the unique 
solution p{t) of X starting at po (that is, the unique curve p{t) satisfying p(t) = X{p{t)) and 
p(0) = Po) is defined and contained in U for all time t G M. 

The problem being local, there are some obvious simplihcations. First, by the classical 
theorem of Darboux, we may assume without loss of generality that (M, H) = (M^”, Hq) where 
Ho is the canonical symplectic structure of and that p* = 0 G Then, we may also 
assume that the one-form ix^o is in fact exact, meaning that X is Hamiltonian: given a 
primitive H of ix^o and letting Jq be the canonical complex structure of the vector 
field can be simply written X = Xh = JqXH, where the gradient is taken with respect to 
the canonical Euclidean structure of M^”'. Therefore 0 is an equilibrium point of Xh if and 
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only if it is a critical point of H, that is VH{0) = 0. Moreover, the Hamiltonian function H 
being defined only modulo a constant, it is not a restriction to impose that H{0) = 0. 

Let (x, y) = (xi,..., x„, yi,..., y„) be symplectic coordinates defined in a neighborhood 
of the origin 0 G so that {x{t),y{t)) = XH{x{t),y{t)) is equivalent to the system 

x{t) = dyH{x{t),y{t)), y{t) = -dxH{x{t),y{t)). 

Since H{0) = 0 but also VLf(O) = 0, the Taylor expansion of H at the origin is of the form 

H{x,y) = H 2 {x,y) + 03 {x,y) 

where H 2 is the quadratic part of H at the origin and where 0^{x, y) contains terms of order 
at least 3 in (x,y). We can now define the linearized Hamiltonian vector field at the origin 
to be the Hamiltonian vector field associated to H 2 '. 

Xh2 = Jo^H2 = Jo A 

where A is the symmetric 2n x 2n matrix (corresponding, up to a factor 2, to the Hessian 
of H at the origin) such that H 2 {x,y) = A{x,y) ■ (x,y). In order to study the stability of 
the equilibrium point, it is useful to first study its linear stability, that is, the stability of the 
origin for the linearized vector field (the latter is obviously equivalent to the boundedness of 
all its solutions). The matrix JqA possesses symmetries which imply, in particular, that if A 
is an eigenvalue then so is its complex conjugate A. It follows that if JqA has an eigenvalue 
with a non zero real part, it also has an eigenvalue with positive real part and in this case one 
can find solutions of the linear system that converges to infinity at an exponential rate: this 
implies linear instability but also instability in the sense of Lyapounov. We will say that the 
equilibrium point is elliptic if the spectrum of the matrix JqA is both purely imaginary and 
simple. This implies linear stability, while linear stability is equivalent to JqA being semi¬ 
simple and its spectrum purely imaginary (but the assumption that the spectrum is simple, 
which is already a non-resonance assumption, will be important for us in the sequel). Note 
that if we only assumed the spectrum to be purely imaginary, then, if the matrix JqA has a 
non-trivial Jordan block, one can find solutions for the linearized vector field converging to 
infinity at a polynomial rate, implying linear instability (but not necessarily instability in the 
sense of Lyapounov). 

So from now on, 0 G is assumed to be an elliptic equilibrium point of the Hamiltonian 
system defined by a smooth function H. Since the spectrum of the matrix JA is invariant 
by complex conjugation, it has necessarily the form {iiai,... ,±ia„} for some vector a = 
(ai,..., a„) G M" with distinct components: this is usually called the frequency vector. By 
a result of linear symplectic algebra (a simple case of a theorem due to Williamson, see 
[AKNOBj i one can find a linear symplectic map which puts the quadratic part into diagonal 
form (this result requires the components of a to be distinct): hence we can assume that H 
is of the form 

n 

H{x,y) = '^aj{x‘j +yj)/2 + 03 {x,y), (1.1) 

i=i 

where our standing assumption from now on is that the Hamiltonian H is real-analytic, hence 
it can be extended as a holomorphic function on some complex neighborhood of the origin. 
Also, we will always assume that the frequency vector a is non-resonant, that is for any 
non-zero A: G Z”, the Euclidean scalar product k ■ a is non-zero. 
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Note that fixing such coordinates imposes a sign on the components of the vector a G M”. 
Given a point (x,y) G let us define I{x,y) G by 

I{x,y) = {h{xi,yi),...,Inixn,yn)), Ij{xj,yj) = (x| + y|)/2, l<j<n 
so that H can be written again as 

H{x,y) = a ■ I{x,y) + 03 (x,y) := hi(I(x,y)) + O^ix.y) 

The linearized vector field, associated to hi{I{x,y)) = a ■ I{x,y), is easily integrated: 
given an initial condition {xo,yo), the corresponding solution {x{t),y{t)) is quasi-periodic. 
More precisely, letting Iq = I{xo,yo) G M”, one obviously has I{x{t),y{t)) = Iq for all time 
f G M and so the set T{Iq) = {{x, y) G | I{x, y) = Iq} is an invariant torus, the dimension 
of which equals the number of strictly positive components of /q, and on which the flow is 
just a flow of translation. The same holds true in fact for an arbitrary Hamiltonian depending 
only on the quantity I{x,y), and such Hamiltonians will be called here integrable. 

A central question in Hamiltonian dynamics is then the following. 

Problem 1. For a Hamiltonian H as in (tni), is the origin stable or unstable? 

By stable we mean Lyapunov stable in the sense that points near the origin remain in a 
neighborhood of the origin. Other notions of stability may also be addressed as we will see 
below. 

1.2 Perturbation of completely integrable systems. 

If H is integrable, the origin is obviously stable. Now in general H is, in a small neighborhood 
of the origin, a small perturbation of the integrable Hamiltonian hi and thus classical tech¬ 
niques from perturbation theory (such as KAM theory, Aubry-Mather theory, Nekhoroshev 
estimates or Arnold diffusion) may be used to tackle the problem. However, this setting of 
singular perturbation theory is quite different from the usual context of a perturbation of an 
integrable Hamiltonian system in action-angle coordinates, that is, a Hamiltonian of the form 
h{I) + sf{9,I), where e is the small parameter and {9,1) G T” x M”. 

A first obvious difference is that for a Hamiltonian H as in (tni, one cannot introduce 
action-angle coordinates on a full neighborhood of the origin: indeed, if we let Ij = Ij{xj,yj), 
then the symplectic polar coordinates 

Xj = ■\/2Ij cos 9j , yj = y^2/j sin 9j, 1 < j < n 

are analytically well-defined only away from the axes Ij = 0. This amounts to the fact that 
for a Hamiltonian integrable in a neighborhood of an elliptic equilibrium point, the foliation 
by invariant tori is singular in the sense that the dimension of each leaf is non-constant (it 
varies from 0 to n), whereas in action-angle coordinates this foliation is regular. 

A second difference lies in the fact that for Hamiltonians of the form h{I) + £f{9,I) the 
perturbation / is usually considered as arbitrary whereas in uni) the perturbation is more 
restricted as it is given by the higher order terms 03 {x,y). 

Finally, a third difference is that, under the assumption that a is non-resonant, a Hamil¬ 
tonian H as in 0 possesses infinitely many integrable approximations /i™', for any integer 
m > 2 (given by the Birkhoff normal form, see below for more details) which are uniquely 
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determined (once the vector a is fixed). This is in sharp contrast with a Hamiltonian of the 
form h{I) +£f{9,I) which does not have, in general, further integrable approximations. 

As we will see below, these differences have the following general effect: in a neighborhood 
of an elliptic equilibrium point, as opposed to a perturbation of an integrable system in action- 
angle coordinates, stability properties are stronger and instability properties are harder to 
exhibit. 

1.3 KAM stability 

Due to the classical KAM (Kolmogorov-Arnold-Moser) theory, one can prove, for any number 
of degrees of freedom and assuming some non-degeneracy assumption (on the higher order 
terms 0^{x,y)), that the elliptic equilibrium point is KAM stable: in any sufficiently small 
neighborhood of the origin, there exist a positive measure set of Lagrangian invariant tori, on 
which the dynamics is conjugated to a linear flow, having the origin as a Lebesgue density 
point. In general, KAM stability does not have direct implications on Lyapounov stability. 

There are however two cases for which one knows that stability holds true for a Hamilto¬ 
nian H as in m- 

The first case is when the quadratic part H 2 is sign-definite, or, equivalently, when the 
components of the vector a G R” have the same sign (and this includes, as a trivial instance, 
the case n = 1). Indeed, the Hamiltonian function has then a strict minimum (or maximum) 
at the origin, and as this function is constant along the flow (it is in particular a Lyapounov 
function) one can construct, using standard arguments, a basis of neighborhoods of the origin 
which are invariant, and the latter property is obviously equivalent to stability. 

The second case is when n = 2 and when the so called Arnold iso-energetic non-degeneracy 
condition is satisfied. Then, KAM stability occurs in every energy level passing sufficiently 
close to the origin, implying Lyapounov stability as the two-dimensional tori disconnect each 
three-dimensional energy level (see for instance [Arnblj and [Mos62] L It is easy to see that 
the Arnold iso-energetic non-degeneracy condition is generic in measure and topology as a 
function of the coefficients of the 04 (x, y) part of the Taylor expansion of H around the origin. 

Related to the results that we will expose in the following sections, let us mention that it 
is sometimes possible to replace the non-degeneracy assumption in the study of stability by 
arithmetic conditions on the frequency vector a of the linear part of the flow at the equilibrium. 
Indeed, in the analytic setting, Herman conjectured the KAM stability (without the Lebesgue 
density requirement) of Diophantine equilibria without any non-degeneracy assumption. In 
f |Her98] l he made the following conjecture (in the slightly different context of symplectic 
maps). 

Conjecture 1 (Herman). Assuming that a. is Diophantine, in any suffieiently small neigh¬ 
borhood of the origin there exists a set of positive Lebesgue measure of Lagrangian invariant 
tori. 

Recall that a G R"" is said to be Diophantine if for some constant 7 > 0 and exponent 
T > n — 1 it holds that |A: • a| > 7 |A:| for all k = {ki,..., kn) G Z” \ {0}, where |/i:|i := 
l^il + • • • + l^nl- We then use the notation a G DC(r, 7 ). 

Herman’s conjecture is true for n = 2, even in the smooth category, as it was proved 
by Riissmann (see for instance |Rusn2| and [FKOQj in the discrete case, for respectively real- 
analytic and smooth maps, and |EFK13j or |EFK15[ Section 7.1] in the continuous case) 
but unknown in general (see [EEKl.Sl IEEK15j for partial results). Note that KAM stability 
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of a Diophantine equilibrium for a Hamiltonian in the case n = 2 does not imply a priori 
Lyapunov stability. 

Observe also that this KAM stability phenomenon without any non-degeneracy condition 
has no counterpart for perturbed integrable system in action-angle coordinates, since any 
integrable system that does not satisfy the so-called Riissmann non-degeneracy condition can 
be simply perturbed so that no invariant torus survives (see |Sev03j i. 

1.4 Arnold’s diffusion conjecture 

Arnold conjectured that apart from these two cases (the case of a sign-definite quadratic 
part, and generically for n = 2), an elliptic equilibrium point is generically unstable. More 
precisely, in [Arn94] one can find the following conjecture. 

Conjecture 2 (Arnold). An elliptic equilibrium point of a generic analytic Hamiltonian sys¬ 
tem is Lyapounov unstable, provided n> 3 and the quadratic part of the Hamiltonian function 
at the equilibrium point is not sign-definite. 

This conjecture is wide open, to such an extent that under our standing assumptions 
(real-analyticity of the Hamiltonian and a non-resonance condition on the frequency vector) 
not a single example is known. 

If the frequency vector is resonant, it is quite trivial to construct an example of unstable 
elliptic equilibrium point (see [MosGOj i. The genericity is, however, still open (see [KMV04] 
for an announcement on some partial results). 

If the Hamiltonian is smooth non-analytic, examples have been constructed by Douady-Le 
Galvez f |DLC83] i for n = 3 and by Douady l |Dou88j i for any n > 3, but here also, genericity 
seems out of reach. 

1.5 Effective stability 

The aim of this paper is to investigate the so called effective stability of an elliptic equilibrium 
point. More precisely, given r sufficiently small and any initial condition (xq, yo) a-t a distance 
at most r from the origin, we are interested in the largest positive time T(r) for which the 
solution {x{t),y{t)), starting at {xo,yo), stays at a distance at most 2r from the origin, for 
all \t\ < T{r). Arnold’s conjecture states that for n > 3, it holds generically that T(r) < oo. 
At the moment there is no other conjectural upper bound on T{r). In this paper, we will be 
interested in lower bound on T{r). Let us first recall some previous results. 

First, without any assumptions, it is easily seen from the equations of motion that T{r) 
is at least of order r~^. Then, given an integer A > 4, with the assumption that H is smooth 
and a is non-resonant up to order K, that is 

kelH, 0< |fc|i < 

the following statement can be proved (see |Bir66j or [Dou88j 1: there exists a symplectic 
transformation 4>^, well-defined in a neighborhood of the origin, such that 

H o ^^{x, y) = a- I{x,y) + h^{I{x, y)) + f^{x, y) (BNF) 

where is a polynomial of degree m = [A/2] (the integer part of A/2) in n variables, with 
vanishing constant and linear terms, and fK is of higher order OK+i{x,y). The polynomial 
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a ■ I{x,y) + h'^{I{x,y)) is usually called the Birkhoff normal form of H of order K. Since 
the term a ■ I{x,y) will be fixed in the sequel we will denote h^{I{x,y)) by BNFk{H). 
The polynomial BNFi^(i/) is uniquely defined, but, in general, this is not the case for the 
coordinate change function (although there is a distinguished choice of a generating 
function for An obvious consequence of (IBNFh is that, in this case, T{r) is at least 

of order at the origin (naturally, the neighborhood in which the effective stability 

holds depends on K and may be very small depending in particular on the arithmetics of 
a). Thus if a is non-resonant and H is of class C°°, T[r) becomes larger near the origin 
than any power of r~^. Observe that if a is non-resonant, one can find a formal symplectic 
transformation and a unique formal series h°° in n variables such that H o ^°°[x,y) = 
h°°{I{x,y)). However, the formal transformation is in general divergent (see [Sie41] l. and 
the convergence problem for the formal series h°° is still an open problem (see |PM03] for 
some results). 

Now with the assumption that the Hamiltonian H is real-analytic, exponentially large 
lower bounds for T{r) have been obtained in two different contexts. 

First, if a is Diophantine, a € DC(t, 7 ), one can prove that T(r) is at least of order 
exp This is obtained by estimating the size of the remainder term in the 

Birkhoff normal form of order K, and then choosing K = K{r) as large as possible in terms 
of r (see |rTDF+89] or |DG96] for slightly better estimates). One should point out here that 


actually for any non-resonant a one can associate a function Aa{r) and prove that T{r) is at 

least of order exp (AQ,(r“^)) (see Section fL 6 l below for the definition of this function A^). In 

1 

the Diophantine case one has Aq,(x) > ( 7 x)'^+i and the classical result is thus recovered. 

Then, in a different direction, assuming only that ol is non-resonant up to order K, for 
some K > but requiring that the quadratic form is positive definite (which implies that 
hi + h?, and then hi + /i™ for any m > 2 , is convex in a neighborhood of the origin), it has 


been proved that T(r) is at least of order exp (r 


K-3 

2n 


this was established independently 


by Niederman l [Nie98| l and Fasso-Guzzo-Benettin ( [FGB98'| 1 and later clarified by Poschel 
f |Pos99j i. The proof is based on the implementation of Nekhoroshev’s estimates l |Nek77] . 
|Nek79| l: observe that in the absence of action-angle coordinates, this implementation is not 
straightforward and it was only conjectured by Nekhoroshev. 

It is a remarkable fact that both exponential stability results under one of the two hy¬ 
pothesis : 1) a is Diophantine or 2) h? is positive definite, can be combined into a double 
exponential stability result if both 1) and 2) hold. This was first done by Giorgilli and Mor- 
bidelli in [MG95j in the context of a quasi-periodic invariant Lagrangian torus. In our context 
of an elliptic equilibrium, the result of |MG95] would amount to double exponential stability 


of a Diophantine equilibrium provided /i^ is positive definite, or more precisely that T(r) is 
at least of order exp ^(exp(( 7 r“^)iT?))^y Even though the condition that h? is positive 
definite is open, it is far from being generic in any sense and recently some efforts have been 
made to improve this result, especially in [Boullb] and |Niel3j . In |Boullbj . using results 
from [Nien7] and [BN12| , it was proved that under a certain condition on the formal Birkhoff 
series /loo, the double exponential stability holds true. This condition, which includes the con¬ 
dition that h? is positive definite as a particular case, was proved to be prevalent (a possible 
generalization of “full measure” in infinite dimensional spaces) in the space of all formal series. 
This result has at least two drawbacks. First, although this condition can be termed generic 
in a measure-theoretical sense, it is far from being generic in a topological sense. Secondly, 
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this condition was only formulated in the space of formal series, and it was unclear whether 
prevalent Hamiltonians have formal Birkhoff series satisfying this condition. This second issue 
was partially solved in |Niel3j : it is proved there that a prevalent Hamiltonian has a formal 
Birkhoff series satisfying a condition close to the one introduced in [Boullb] , yielding a result 
which is only intermediate between exponential and double exponential stability. 

The aim of this paper is to improve those results by establishing that generically, and in 
a strong sense, the double exponential stability holds true. 

1.6 Main results 

We start by some reminders and notations that will be useful in our statements. Let H be 
a real analytic Hamiltonian on having an elliptic equilibrium point at the origin with a 
non-resonant frequency vector a, that is H is as in (HI]). 


• For vectors in || . || denotes the norm defined as 


z := max a 

J P ? — ('^1; ■ ■ ■ ? • • 

to 

(1.2) 

l<j<n 



and for vectors in C”', . 

denotes the usual Euclidean norm 



I|/|| := 

V|/l|2 + --- + |/„|2, /= (Ii,. ..,/„,). 


(1.3) 


It will be more convenient to use these different norms for vectors in or in C"', and 
we hope that this abuse of notations will not confuse the reader. 


• We suppose that the radius of convergence of H is strictly larger than some R > 0 and 
let ||H||/{ be the sup norm of H in the open complex ball in centered at the origin 
of radius R that we denote by 

Br := {z e I ||z|| < R}. (1.4) 

We also define the real ball Br := Br n 

• We denote by P{n,m) the set of polynomials of degree m in n variables. We let 
P 2 {n,m) C P{n,m) be the subspace of polynomials with a vanishing affine part, and 
P‘i{n, m) C H(n, m) the subset of polynomials that have a vanishing affine and quadratic 
part. 

• We denote by Hm G P3{2n,m) the part of the power expansion of H that contains the 
terms of degree between 3 and m included. 

• Having fixed the number of degrees of freedom n, in all the sequel, we let 

Ko = Ko{n) := +4, niQ = mo{n) := [Ko{n)/2]. 

• The vector a is supposed to be non-resonant: this means that for any integer K > 1, 

= max{|A: • a\~^ \ k G Z”, 0 < |A:|i = \ki\ + ■ ■ ■ \kn\ < K} < -|-oo. (1.5) 

We define, as in |Boul2] . the function 

^a{x) = SUp{iL > 1 I K'I/a{K) < x}. 
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Observe that if a G DC(t, 7 ), then '^aiK) < 7 and hence 

> (7x)l+^ (1.6) 

• Recall that for H as in o, there exists for every integer RT > 4 a real analytic 
symplectic transformation < 1 >^ defined in the neighborhood of the origin such that 

H o^^{x,y) = a ■ I{x,y) + h'^{I{x,y)) + f^{x,y) 

where is a polynomial of degree m = [K/2] (the integer part of K/2) in n variables, 
with vanishing constant and linear terms, and /k is of higher order OK+i{x,y). We 
denoted h™ by BNFk{H). By uniqueness of the Birkhoff normal form we have for 
K = 2m > 4, a well defined map 

BNFi^ : Ps{2n,K) —> P 2 {n,m) 

Hk = BNFK(^i^) = BNFk{H). 

Our main result is the following. 

Theorem A. Let H be a real analytic Hamiltonian on having an elliptic equilibrium 
point at the origin with a non-resonant frequency vector a. There exists an open and dense 
set of full Lebesgue measure A/’n(a) G P 3 ( 2 n,R'o) such that if Hkq G J\fn{a), then there exists 
r*,c,c',c'' > 0 that depend only on n,R,\\H\\ii,a and Hkq such that if r <r*, then 

T{r) > exp (cr“^exp (c'Aq (c'V“^))) . 

If a £ DC(t, 7 ), there exists an open and dense set of full Lebesgue measure A/’n(«) £ 
P 3 ( 2 n,R'o) such that if Hkq G Mn{oi), then there exists r* and C that depend only on 
n,R, and Hkq such that if r < r*, then 

T{r) > exp ^exp . 

Observe that since d and d' will not depend on a (see ()2.8l) i. it follows from (ll. 6 p that 
the constant C that appears under the double exponential in the Diophantine case is actually 
of the form C = 7 "^O' where C does not depend on a. Theorem lAl improves all previous 
results contained in [MG95j . [Boullb] and |Niel3j . In the course of its proof, we will also 
have to extend the results on exponential stability contained in |Nie98| . [FGB98| and [Pds99j . 

Remark 1.1. Observe that even though AQ,(r“^) goes to infinity as r goes to zero, the speed of 
convergence can be arbitrarily slow but the statement implies that T{r) is always at least of 
order exp(cr“^). From the proof of the theorem, one can easily obtain the following statement: 
fixing /c G N*, k > 2, and allowing the constants and to depend also on k, one has 

T(r) > exp (^Ckr~^ exp (c'A^ (c'V“^))^ 

which is always at least of order exp(cfcr“^). As a matter of fact, the weaker estimate 

T{r) > exp (^Ckr~^^ 

can be obtained if one only assumes a to be non-resonant up to a sufficiently high order 
depending on k and n. 
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Remark 1.2. The Diophantine condition a G DC(t, 7 ) is sometimes called an asymptotic 
Diophantine condition. A strictly weaker condition, called uniform Diophantine condition, 
requires the existence of an increasing sequence Kj G N, Kj —)■ 00 , such that |A: • a| > 
for every A: G Z” \ {0} with \k\i < Kj. This gives ^a{Kj) < 'jKJ and Theorem lAl would then 
imply that there exists a sequence rj —)• 0 such that 

T{rj) > exp I exp 



The notion of stably steep polynomials, which can be implicitly found in the work of 
Nekhoroshev 1 |Nek73| l. will be important in the proof of Theorem lAl 

Definition 1 (Stably steep polynomials). A polynomial Pq G P 2 {n,m) is called stably steep if 
there exist a neighborhood V of Pq in P 2 {n,m) and positive constants C,6 such that for any 
integer I G [l,n — 1], any P G V and any vector subspace A C M” 0 / dimension I, letting Pa 
be the restriction of P to A, the ineguality 


max 

o<v<^ 


min 

||a:||=r;, xGA 


||VPA(a;)|l>Cr"' 


holds true for all 0 < < 6, where 11 .11 is the usual Euclidean norm defined in m- 

The set of stably steep polynomials in P 2 {n,m) will be denoted by SS{n,m). 

Theorem El will clearly follow from the combination of the following two statements. 
Theorems iBl and O with the set A/’n(a) being defined as Mn{(y) := BNF)(-J(55(n, mo)). 

Our first statement is that the set of Hamiltonians with stably steep BNF of order Kq 
have doubly exponentially stable equilibria. 

Theorem B. Let H be a real analytic Hamiltonian on having an elliptic equilibrium 
point at the origin with a non-resonant frequency vector a. If 


BNFxo(^) = h^° G SS{n,mo) 
then the conclusions of Theorem \A\ hold. 

The second statement shows that the condition BNFKq{H) = BNFKoiHKo) £ SS{n,mo) 
is generic in a strong sense. 

Theorem C. For any non-resonant a G M”, the complement of BNFj}{SS{n, mo)) in 
P3(2n,iFo) is contained in a semi-algebraic subset of positive codimension. In particular, 
BNFj}{SS{n,mo)) is a dense open subset of Po{‘2,n,Ko) of full Lebesgue measure. 

Proof of Theorem HI Putting together Theorem [B] and O immediately yields Theorem H] if 
we take A/’n(a) = BNF~^^^{SS{n,mo)). □ 

To prove Theorem o we will show that the complement of SS{n,mo) in P 2 {n,mo) is 
contained in a semi-algebraic subset of codimension at least one. This will be done in Sections 
12 . 1 [ 12.21 and Appendix HI 

Theorem [B] will follow (see Section El from a version of the Nekhoroshev exponential 
stability result adapted to our singular perturbation setting, that we now present and that 
will be proven in Section [3l 
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• For vectors in C"", it will be convenient to also use the sup norm | . | defined as 

|/| :=max{|/i|,...,|4|}, / = (/i,..., 4). (1.7) 

This norm allows an easier comparison between I{z) G C” and G indeed, we 

have \I{z)\ < ||z|p/2 and the equality holds true if z G 

• Given r > 0, we define the domain to be the open ball centered at the origin in C”' 
of radius r^/2 with respect to the norm | . |: 

Vr := {/ G C" I |/| < r‘^/2} 

and we let Dr := Dr f) M”. This choice is motivated by the fact that if / : z G i—)• 
I{z) G C"", then I{Br) C Vr and I{Br) = Hr H M”, where Br and Br have been defined 
in p.4|) . 

• We define || . ||r to be the sup norm for functions defined on Br or on Vr- Extending the 
norm || . || initially defined for vectors in C"" and (respectively in (II.2p and in (II.3M 
to tensors in C”' and we extend the sup norm || . ||r for tensor-valued functions 
defined on Br or on Vr- The same notation || . ||r will be used also for the real domains 
Br and Dr'- this will not cause confusion as it will be clear from the context if it is the 
complex or the real domains that are considered. 

• We consider a Hamiltonian H of the form 

H{z) = h{I{z)) + fiz), h:Vr^C, f:Br^C {*) 

which is real analytic and such that 

\\Vh\\r<E, \\V‘^h\\r<F, ||X/||r<e (1.8) 

where Xf is the Hamiltonian vector field associated to / 

• The integrable Hamiltonian h is supposed to be steep on the domain Dr, as defined 
below. 

Definition 2. A differentiable function h : Dr M is steep if there exist positive 
constants C,6,pi, for any integer I G [l,n — 1], and k such that for all I G Dr, we 
have ||V/i(/)|| > k and, for all integer I G [l,re — 1], for all vector space A G ffi” of 
dimension I, letting A = / + A the associated affine subspace passing through I and h\ 
the restriction of h to X, the inequality 

max min \\Vh^{r) - Vhx{I)\\ > 

0<V<i ||/'-/||=r?, J'GAnDr 

holds true for all 0 < f < S. We say that h is {r,K,C,6, {pi)i=i^...^n-i)steep and, if all 
the Pi = p, we say that h is (r, k, C, 6, p)-steep. 

Let us point out that the definition of steepness that we use is not exactly the one given by 
Nekhoroshev but it is obviously equivalent to it (see [Nek73| or |Nek77j l. Indeed, Nekhoroshev 
only requires steepness for subspaces A which are orthogonal to V/i(/), in which case Xhx{I) = 
0; for subspaces A such that Vh\{I) / 0, the inequality in Definition [5] is clearly satisfied 
(and one may even set p/ = 0 in this case). 
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Theorem D. Let H{z) = h{I{z)) + f{z) be as in (jlj) satisfying (II.Sp . such that h is 
{r,K,C,6,{pi)i=i^,,,^n-i)-steep. Then there exist r*,c,c' > 0, which depend only on n, E, 
F, K, C and pi for 1 < I < n — 1 such that if 

r<h\ re < (1.9) 


where 


a:=l+pi +piP 2 H- VpiP2 ■ ■ -Pn-i, 


then for any solution z{t) of the Hamiltonian flow (jlj) with z(0) = zq € i?j ,/2 we have 


\I{z{t)) — I{zo)\ < c' {r£)2na ^ |i| < exp ^(re) 2na^ . 


Remark 1.3. We will prove in fact a slightly more general and more precise statement (but 
whose formulation is also more cumbersome): there exist positive constants ci, C2, C3, C4, C5, ce 
and C 7 , which depend only on n, E, F and on the steepness constants k, C,pi, for 1 < I < n — 1, 
such that for any solution z{t) with z{0) = zq ^ Br/ 2 ^ if 


re < min 



,C3r 


Ana 



( 1 . 10 ) 


where a is as above and 


a' ;= 1 + p 2 + P 2 P 3 H-h P 2 P 3 ■ ■ ■ Pn-i, 

then 

\I{z{t)) — I{zo)\ < c^{re)^^, |t| < C6r“^(re)“^ exp ^C 7 r“^(re)“^^ . 

This statement obviously implies the statement of Theorem [Dl Let us also add that using 
this more precise statement, one can easily obtain a more precise statement in Theorem K1 

1.7 Comments, open questions and prospects 

It is natural to ask whether our main result, Theorem IXl can be improved, and so we can ask 
the following two questions. 

Question 1. Does Theorem\^ remains true without assuming BNFxo(17i<'o) ^ SS{n,mo)? 

Question 2. Is the estimate on the time T[r) in Theorem\^ essentially optimal? 

A main difficulty in these questions is related to the fact that the construction of an 
unstable elliptic equilibrium point in the analytic category is a wide open problem as we 
emphasized in the Introduction. Concerning the second question, let us just mention that it 
may be possible to give an answer in the Gevrey category (a regularity which is intermediate 
between smooth and analytic). Indeed, on the one hand, one should expect that the statement 
of Theorem [A] holds true for Gevrey Hamiltonians, with only different constants. On the 
other hand, using the methods in |MS02j . it might be possible to construct an unstable 
elliptic equilibrium point in the Gevrey category, with a time of stability which is a double 
exponential (the fact that one can construct a Lyapunov unstable elliptic equilibrium point in 
the Gevrey category follows directly from |Dou88] , but the real difficulty is to get an estimate 
on the time of instability). 
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Then it is also natural to ask whether our result holds true for a quasi-periodic invari¬ 
ant Lagrangian torus, or more generally, for a quasi-periodic normally elliptic and reducible 
invariant torus (which includes both elliptic equilibrium points and quasi-periodic invariant 
Lagrangian tori as particular cases). This general case is described by a Hamiltonian of the 
form 

H{e,J,x,y) = j3 ■ J + a- I{x,y) + F{e,I,x,y) 

where {9, J) £ T™' x M™' are action-angle coordinates, {x, y) symplectic coordinates around 
the origin in and F is at least of order 2 in / and 3 in {x,y). The set {{J,x,y) = 0} = 
{(J, x,y) \ J = 0, I{x, y) = 0} is a normally elliptic torus of dimension n in a n -|- m degrees 
of freedom Hamiltonian, and the question is as follows. 

Question 3. Assuming that the vector (/3, a) G is Diophantine and H is real-analytic, 

does Theorem extends to this setting in the following sense: i/(J(0),/(x(0), y(0))) is at a 
distance r of zero in with r sufficiently small, is it true that {J{t),I{x{t),y{t))) stays 

at a distance 2r from 0 for a time T{r) which is doubly exponentially large with respect to 
7 .-i/('r+i) fyjhere r is the exponent of the Diophantine condition on the vector {fd,a))? 

In a subsequent paper 1 [BFN15] 1. we will answer positively the above question in the 
case of an invariant Lagrangian Diophantine torus that is of particular interest in the study 
of perturbed integrable systems. Indeed, by KAM theory, it is well-known that invariant 
Lagrangian Diophantine tori appear for arbitrary small perturbations of generic integrable 
Hamiltonian systems in action-angle coordinates. Furthermore, these tori are not isolated and 
appear as a family parametrized by some Cantor set of positive Lebesgue measure (tending 
to full measure as the size of the perturbation goes to zero). The goal of |BFN15j is to prove 
that under an additional generic assumption on the integrable Hamiltonian, most of the KAM 
tori are doubly exponentially stable. 

2 Genericity of steepness and Birkhoff normal forms 

The aim of this section is to give a proof of Theorem O and of the fact that Theorem [Pl 
implies Theorem [Bl 

2.1 Genericity of steepness 

In Appendix lAl we will prove a general result on genericity of stably steep polynomials. 

Theorem 2.1. The complement of SS{n,mo) in P 2 {n,mo) is contained in a semi-algebraic 
subset T(re, mo) of codimension at least one. 

Theorem 12.II has an immediate consequence on the genericity of steep functions as will be 
shown in the following Theorem 12.21 

Given p G N, p > 3 and p > 0, let C^{Dp) be the set of functions p times continuously 
differentiable on Dp, and let 

ll^'^^llp < oo 

where || . ||p is the sup norm on Dp of the tensor-valued function and where we recall 

that by definition, Dp is the (real) open ball of radius p ^/2 with respect to the sup norm | . |. 
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Given h G C^{Dp), we denote by Tp-ih{I) G P{n,p — 1) the Taylor expansion of h of order 
p — 1 at / G Dp (or the p — 1-jet at I). We have the following statement (that will be used 
later with the value p = mo -|- 1). 

Theorem 2.2. Let h G C^{Dp) be such that ||V/i(0)|| := w and Pp-i := Tp-ih{0) — Tih(O) — 
Toh{0) G SS{n,p— 1). Then, there exists positive numbers p*,S* and C that depend only on 
w, Pp-i, ||V^/i||p^p and n such that h is {p, K,C,d,p — 2)-steep, with 

p := mm{p/2, p*}, k ■.= w j2, 5 := min{p^/4, (5*}. 

Proof of Theorem \2.iA Let M := ||V^h||p^p. Observe first that if p^ < w/M then the condi¬ 
tion ||V/i(/)|| > K = wI2 is satisfied for any 1 G Dp. 

Fix an arbitrary 1 G Dp, and define TLj = Tp-ih{I) — Tih{I) — Toh{I) G P 2 {n,p — 1). 
Since LLq = Tp-ih{0) — Tih{0) — Toh{I) = Pp-i is stably steep, we have the existence of p 
that depends on M, Pp-i, and n such that p < p, Tii is sufficiently close to Pp-i so that 
for all integer I G [1, n — 1], for all vector subspace A C of dimension I, letting be the 
restriction of "H/ to A, the inequality 

max min ||V^/,a(x)|| > 

holds true for all 0 < ^ < (5o, where (5o and Cq are the steepness constant related to Pp-i. 
Now, we get by the Taylor formula (applied to V/i at the order p — 1) that 

\\Vh{I + x)- Vh{I) - Vni{x)\\ < M{p - l)!||x||P-i 

provided I -|- a; G Dp, which is satisfied if p < pl2 and |x| < ||x|| < p^/4. So for ||x|| < f, < 5, 
with 5 := min{Co(2M(p — , p^/A}, we have 

iivh(/+x) - v/i(/) - vHi{x)\\ < {Come-"^ 

and then, letting A = I -|- A, 

||V/ia(/ + x)- Vhxil) - VnpA{x)\\ < {Co/2)e-^. 

From this we eventually obtain 

max min \\VhxiI + x) - V/ia(/)|| > {Co/2)£f’-‘^ 

0<rj<^ \\x\\=r],x£A 

and letting I' = / -|- x, C := C'o/2, 6* := Co{2M{p — 1)!)“^ and p* := min{/i, vj/M}, the 
steepness of / is thus established with the constants given in the statement. □ 

2.2 Generic steepness of the BNF. 

The proof of Theorem O will be an easy consequence of Theorem 12.11 and the following two 
lemmas on the map BNF a:. 

Lemma 2.3. The map BNFa' is algebraic. 

Proof. This follows by construction of the Birkhoff normal form, and we refer to |PM03j for 
more details. □ 
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Now given a polynomial Q = Q2 + ■ ■ ■ + Qm £ P2{n, m), where each Qj is homogeneous 
of degree j, it can be identified to a polynomial Q G P‘i{2n,K) by setting 0(0 •= 

For iF > 4, we can define a map by 

: P2{n,m) — > P2{n,m) 

Q ^b1^Fk{Hk + Q)- 

Lemma 2.4. The map preserves Lebesgue measure. 

Proof. This also follows by construction of the Birkhoff normal form. More precisely, it 
can be shown that decomposing the map as where Ff is the 

component with respect to homogeneous polynomials of degree j, then we have F^{Q) = 
Q 2 + BNF\H 4 ) = Q 2 + h\ and for 3<j<K,we have Ff{Q) = Qj+Pf {H 2 j,Q 2 , ■■■, Qj-i) 
where pP is an algebraic map (see (Nie m, where this property has already been used). 
This expression clearly implies that is smooth with Jacobian one, therefore it preserves 
Lebesgue measure. □ 

Proof of TheoremK\ Our aim is to show that the complement of (S'5'(n, mo)) in 

P3(2n,iFo) is contained in a semi-algebraic subset of positive codimension. Since the inverse 
image of a semi-algebraic subset by an algebraic map is semi-algebraic, from Theorem 12.II and 
Lemma 12.31 it follows that the complement of BNF)^^(S'S'(n, mo)) in P^[2n, Kq) is contained 
in a semi-algebraic subset. It remains to prove that this set has positive codimension, or equiv¬ 
alently, zero Lebesgue measure in P3{2n, Kq). By Lemmafor any Hkq £ -^ 3 ( 2 ^, JCo)) the 
Lebesgue measure in P 2 {n,mo) of the set 

{Q G F 2 (n,mo) \Hko + Q^ BNF]^;(5S(n, mo))} 

is zero. By Fubini-Tonelli theorem, this implies that the complement of BFfFj}^{SS{n,mo)) 
in F3(2n, Kq) has zero Lebesgue measure, and this concludes the proof. □ 


2.3 Birkhoff normal forms with estimates 


For a real analytic Hamiltonian with an elliptic equilibrium point, as in m, it is known 
that the estimates on the Birkhoff normal form are given by the arithmetic properties of a 
and the analytic norm of H. We summarize in the following Proposition 12.51 the estimates 
on the BNF that will be useful for us in the sequel. The proof of Proposition 12.51 is relatively 
standard, we include it in Appendix [B] following |DG96j . 

Here it will be more convenient to perform a linear change of complex canonical coordinates 
■2 = <5(0, where S : —)• is dehned by 






^ (0 '^^n+j ) • 


It is easy to check that this linear transformation S and its inverse S~^ have unit norm 
(with respect to the norm || . || defined in (II.2M . hence H and Ff o S have the same radius of 
convergence around the origin and ||Ff o Abusing notations, we will still write 

H instead of Ff o S to denote the Hamiltonian in these new coordinates. Observe that 


if2(0 = hiim) = a • /(O = f ^ ajCjCn+j 

t=i 
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where 


I{C) = Ij{0=iCjCn+j, l<j<n. 

Recall the definition of 'I'q given in Section 11.61 and define also for any integer j > 3, 

3 

i=3 


For RT > 1, define 

PK ■■= {548ncdK^{K))-^, (2.1) 

where the positive constants c and d depend only on n, R and ||^f||_R and are defined in ()B.3p . 

Proposition 2.5. Let H be as in dEU) with a as in (lESD, and fix integers p > 2, K > 
2p and 0 < q < K — A. There exist constants b{p) and b{q) that depend respectively on 
p, n, R, ||i?||_R, and on q, n, R, ||^f ||_r, such that if we assume 

0 < p < p^/e, (2.2) 

then there exists a real-analytic symplectic transformation = Id + O(C^) defined on Bpj^ 
such that 

H o ch^(0 = a • 1(0 + h^{I{0) + f^iO, 
with f^ = 0 (^^+^) and the following estimates hold 

WV^h^Wpp = max ||V^'/i'"|L < b{p) (2.3) 

2<j<p 

\\Vf^\\p<b{q)p'^e-^. (2.4) 


2.4 Prom Nekhoroshev stability to double exponential stability 

In this section we prove that Theorem [D] implies Theorem [Bj As a corollary of Proposition 
12.51 and Theorem 12.21 we get the following 

Proposition 2.6. Let H be as in (ini) with a as in CSI), and such that 

k^Hk,) = G R5(n,mo). 

There exists C > 0 and K* > 4 that depend only on n, R, ||^I||_r, hmQ, lla|| and such 

that if 

K > AT*, ^ < p < Pk/^, 

then there exists a real-analytic symplectic transformation = Id + 0(^^) defined on Bpj^ 
such that 

H o $^(0 = a • 1(0 + h-(/(0) + f^{^) ■■= h{L{0) + f^{0 (2-5) 

with f^ = 0 (^^+^) and 




\mo+l,p 


max IL < 6(mo + 1) 
2<j'<mo+l 


and such that h is {p/2, 


\\Vf^\\p<b{q)p^e-^, 0<q<K 

a||/2, C, p^/4, mo — l)-sfeep. 


-4, 


( 2 . 6 ) 

(2.7) 
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Proof of Provosition \2.(k For K > K* > 2(mo + l) apply Proposition [23] withp = rriQ + l and 
q < K — A and get (12.51) with estimates (|2.6p and (|2.7I) . We want to apply Theorem 12.21 with 
p = mo + 1 and vo = ||al|. Observe first that Tmoh{0) — Tih{0) — To/i(0) = /i"*° € SS{n,mo). 
Then observe also that V^/i = V^/i™ and that by (12.61) . we have the bound 

||V^/i||mo+i,p = ||V^h™||mo+i,p < + 1) 

which is independent of p, hence the constants C, p* and <5* in the statement of Theorem 12.21 
do not depend on p, and choosing K* sufficiently large, pK and then p become sufficiently 
small so that pl2 < p* and p'^/A < 6* therefore h is {p/2, ||a| |/2, C, p^/4, mo — l)-steep. □ 

We now use Proposition 12.61 and Theorem iDl to give the 

Proof of Theorem\^ Let H be as in dLH) with a as in (11.51) and 

BNFko{Hko) = G SS{n,mo). 


For r > 0 we define 

K = Aa{{l6AAencdr)~^) 

so that pkI& > 3r, and observe that K > K* is satisfied (with K* given by Proposition 12.61) 
provided r < r* for some sufficiently small r* > 0. Hence we can apply the latter proposition 
with our choice of K and with p = 3r. 

Next we want to apply Theorem iDl to ()2.5I) . First observe that Theorem El is stated and 
proved in the z variables whereas the estimate of Proposition 12.61 are given in the ^ variables: 
however since z = S{f) with S and 5“^ of unit norm, Theorem iDl also holds true, with the 
same estimates, if one uses the ^ variables. 

From Proposition 12.61 and our choice of p, the function h is (3r/2, ||a| |/2, C, 9r^/4, mo — 1)- 
steep and (jl.Sp is satisfied with 

E:=3/2||al|, F:=b{2), e := b{q)p‘ie-^ 

for some 0 < q < K — A yet to be chosen. Up to taking r* smaller one easily checks that dLlI) 
(with r replaced by p = 3r) is satisfied provided we choose q = Ana — 1. Thus Theorem iDl 
can be applied and we obtain the following statement: given an arbitrary solution ^(t) of the 
system associated to H o in (12. 5p . if ||f(0)|| = 1|5(0)|1 < p/2 = 3r/2, then 

\Iiz{t)) - /(z( 0 ))| < c'ipe)^, |t| < exp (^{pe)~^'^ . 

For r sufficiently small, this implies in particular that ||C(t)|| = ||5(i)|| < 7r/4 for times 

\t\ < exp ^(3re)“^^ . 

Recalling the definition of e and with our choices of q and K, the previous estimate implies 
that ||f(t)|| = ||5(t)|| < 7r/4 for times 

\t\ < ex.-p{cr~‘^ exp{c'Aa{c''r~^))) 

with 

c := 3“^6(4na — l)~^, d := (2na)“^, d' := IGAAencd. (2.8) 
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To conclude, observe that H is related to (12.Sp by a symplectic transformation = 
Id + which can be made close enough to the identity (as well as its inverse) by taking 

r small enough. Hence, given any solution ^(t) of the system associated to H with ||^(f)|| = 
lk(i)|| < r, the corresponding solution ({t) of (12.511 satisfy ||C(0)|| = ||5(0)|| < 3r/2 for r 
small enough, and therefore ||^(t)|l = l|^(t)|| < 7r/4:, and also ||^(t)|| = || 2 (f)l| < 2 r, for times 

\t\ < ex.p{cr~'^ exp{c'Aa{c"r~^))). 

We eventually arrives at the estimate 

T{r) > exp(cr“^ exp(c'AQ,(c'V“^))) 

and this concludes the proof of the lower bound on T{r) in the general case. The estimate in 
the Diophantine case follows from the general case and from (II. 6 p . □ 

3 Nekhoroshev exponential stability for an elliptic equilibrium 

The goal of this section is to prove Theorem [Dl following the method introduced in |BN12j 
and [Boullaj . This method, which uses only periodic approximations and compositions of 
periodic averagings, has the advantage of being directly applicable in a neighborhood of an 
elliptic equilibrium point where action-angle coordinates cannot be used. 

Since the proof contains some technical statements, we first give in the next Section [3.11 a 
long and complete heuristic description of the method that would hopefully make the reading 
of the proof easier. We emphasize that Section [3. II is included only for the convenience of the 
reader and does not interfere with the proof strictly speaking. 

3.1 Heuristic description and plan of the proof 

Given an arbitrary initial condition zq and the associated solution z{t) (that is 2 :( 0 ) = zq) of 
the Hamiltonian H = h + f, our goal is to prove that the variation of the action I{z{t)) — I (zq) 
remains small (as a small power of e) for an interval of time which is exponentially large with 
respect to the inverse of (some power of) e, £ being the size of the perturbation /. The proof 
is based on an algorithm that, for 0 < j < n — 1 , reduces to a space of dimension n — j — 1 the 
directions in which a fast drift (before an exponentially long interval of time) may be possible 
at each step j, and that stops therefore after at most j = n — 1 steps. We now describe the 
heuristics of this algorithm, which depends on a positive parameter Q > 1 and an integer 
parameter m > 1. 

For the step j = 0, we write H = Hq and given the parameter Q > 1, we use Dirichlet’s 
box principle to approximate the unperturbed frequency vq = Vh{I{zo)) by a periodic vector 
ujQ, that is, a vector which is a real multiple of an integer vector (this corresponds to a 
frequency vector which is maximally resonant, as the set of integer vectors k orthogonal to 
ujQ forms a sub-module of maximal rank n — 1). Letting Tq be the period of wq, which is the 
smallest positive number t such that tujQ € Z"", the parameter Q controls the approximation 
as follows: 

Ibo -^^oll = So < {ToQ)~\ Iboir^ ^ ^0 < 

Then, on some small neighborhood Vq oi zq, given the integer parameter m > 1 and assuming 
certain compatibility conditions between m, sq, Tq and e, it is possible to construct a resonant 
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normal form (with respect to a;o) up to a remainder which is exponentially small in m: more 
precisely, by a symplectic transformation $0 which is close to the identity, the transformed 
Hamiltonian Hq o $g can be written as a perturbation of h, but this time the perturbation 
splits into two parts: a resonant part, which is still of order e but has the additional property 
that its Hamiltonian flow commutes with the linear flow of frequency wg, and a non-resonant 
part which is of order As a consequence, we have the following partial stability result 

for the solution of Hq o <bg starting at ‘hg ^(zg): the variation of the action variables in the 
(one-dimensional) direction given by loq is small for an exponentially long interval of time 
of order 2”^, unless the solution escapes from the domain of the normal form Vg before. In 
other words, we excluded at this step j = 0 the direction ujq from the directions along which 
a drift in the actions may appear before an exponentially long interval of time. Since ^g 
is symplectic and close to the identity, the same holds true for the solution z{t) = 2 g(t) of 
H = Hq, and we arrive at the following dichotomy: either the action variables have also small 
variation in the direction transverse to wg for an exponentially long interval of time, or not. 

In the first case, for an exponentially long interval of time with respect to m, the variation 
of the action variables is small: the stability condition is satisfied and the algorithm stops. 
Once the algorithm stops, one can determine the parameters Q and m in order to fulfill the 
compatibility conditions which essentially read as follows: 

So < 1, mTQE < sg, rriToSQ < 1. 

Since sg ~ {TqQ)~^ and Tg < (as ||V/i(/(^;g))|| = ||ug|| is of order one), these conditions 

are satisfied if we choose m ^ Q ^ and e < 1, and thus we obtain a result of exponential 
stability with respect to m ^ If /i is convex or quasi-convex it is simple to see, due 

to energy conservation, that this first case is automatic and it is not possible that the action 
variables drifts transversely to cug, hence exponential stability is proved with one step of the 
algorithm. But in general, the second case is possible and more work is needed to further 
reduce the drifting possibilities of the actions going from the step j = 0 to the step j = 1 of 
the algorithm. 

In the second case, setting Hq = Hq o ^q and denoting by Zq (t) the associated solution, 
we can find a positive time Iq, which is shorter than 2™, such that the I{zQ{iQ)) — /( 2 :()'( 0 )) 
has a small drift of order sg in the direction orthogonal to ujq. Letting Hg be the projection 
onto the orthogonal of uq, we can define a curve 

7o(() = /(Zo+O)) + n„(/(4{()) - 7(4(0))) 

which takes values in an affine subspace of dimension n — 1. One can then exploit the steepness 
property to find a time tg < tQ for which the vector V/i( 7 g(tg)) is linearly independent from 
loq in a quantitative way: 

no(V/7(7o(io))) > 

where Pn-i is the steepness index in dimension n — 1. Using Dirichlet’s box principle again 
with the same parameter Q, we can approximate the vector ng(V/i( 7 g(tg))) = vi by another 
periodic vector uji: 

11^1 - a;i|| = Si < {T,Q)-\ llniir' <T,< llmll-iQ"-' < 

First observe that since vi is orthogonal to ujq, ui is almost orthogonal to ojq and in particular it 
is linearly independent from wg. Then since || 7 g(tg)—/( 2 ;(J'(tg))|| = ||nf|-(/(Z(j'(tg))—/(2:(j'(0)))|| 
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is small (as I{zq {t) — I{zq (0)) has only small variation in the direction given by wq), we also 
have 

||no(Vh(/(4(to)))) -a;i|| ~ si < {T^Q)-\ Ti < 

Set zi = ^(J~(to) and Hi = . On a small neighborhood Vi of zi (small enough so that Vi is 

still included in Vq) we can then construct, as in the first step, a resonant normal form with 
respect to cui up to an exponentially small remainder with respect to m. Unlike the step j = 0 
in which the perturbation was arbitrary, here the perturbation is given by the non-resonant 
part with respect to wQ) and this explains why it is sufficient to have an approximation of 
no(Vh(I( 2 i))) and not of the full vector Vh{I{zi)). Moreover, a careful construction of 
the new normalizing transformation <f>i shows that the resonant part of Hi o <f>i, whose flow 
commutes with the linear flow of frequency ui , also commutes with the linear flow of frequency 
ojQ. Exactly like in the first step, we arrive at a dichotomy which determines whether the 
algorithm stops (the variation of the action of the solution of and then of Hi = Hqo^q 

and Hq = H, is small and the theorem is proved) or moves to the next step with the gain 
that now both the directions ujq and oji, that are linearly independent, are excluded from 
the directions along which a drift in the actions may appear before an exponentially long 
interval of time. Note that if the algorithm stops, the parameters Q and m have to be chosen 
according to 

Si < 1, mTi£ < Si, mTiSi <1, 0 < i < 1. 

Observe that 

Ti < ~ {TuQY-^Q^-^ < Q-(i+Pn-i)-i 

1 

and hence the compatibility conditions are satished if we choose m j'Nj £ e ^ 1, 

with oi = 1 + Pn-i, and thus we obtain a result of exponential stability with respect to 

1 

m ^ £ . In particular, as 2nai > 2n, this stability result also holds true if the algorithm 

stopped at step j = 0. 

We have briefly explained how to pass from the step j = 0 to the step j = 1, and how 
the parameters Q and m are chosen if the algorithm stops at step j = 0 or j = 1. But of 
course, for any 0 < j < n — 2, one proceeds exactly the same way to go from step j to step 
j + 1. The fact that the algorithm stops after n steps (if, of course, it didn’t stop before) 
is clear since then n linearly independent directions are excluded from the directions along 
which a drift in the actions may appear before an exponentially long interval of time. More 
formally, in the case j = n — 1, the resonant part in the normal form HY_i = H^-i o 
consists of a Hamiltonian whose flow commutes with n linearly independent linear flows with 
frequency wq, ..., Wn-i; it is plain to see that such a Hamiltonian is integrable, so HY_i consist 
of an exponentially small perturbation of some integrable Hamiltonian: the first case of the 
dichotomy thus holds, the algorithm stops and the theorem is proved. At each step j, the 
compatibility conditions are given by 



Si < 1, 

mTi£ < Si, mTiSi 

<1, 

0 < i 

VI 



and using the fact that 








To < Q^-\ 

Tj < {Tj_iQr--^Q 

n—1 

1 < j 

< n — 

1, 


we can choose. 

if the algorithm stops at step j, m ^ 

Q ~ 

1 

^ 2naj 

and e 

<1: 

, with 

Oo = 

= 1, ai = l+pn- 

— 1 , dj — 1 + Pn—j 

+ ••• 

“t“ Pn—j 

■■■Pn- 

-1, 

J >2, 
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leading to a result of exponential stability with exponent 2naj. We have 

an-l = l+ Pl+ PlP2 H-h PlP2 ■ ■ ■ Pn-l = a, 

and since 2na„_i > 2 nan -2 > • > 2nao, independently of the step j at which the algorithm 

stops (and hence independently of the choice of Q and m), we obtain a result of exponential 
stability with exponent 2na. 

Let us now describe the plan of the proof. In US.21 our aim will be to obtain a suitable 
normal form (Proposition [3]T]) for an abstract Hamiltonian Hj, where 0 < j < re — 1, which, 
as we explained in the heuristic description above, will be later related to our original Hamil¬ 
tonian H described in Q in the following way: Hq = H and for j > 1, Hj = Hj-i o = 
Hq o $g o • • • o where <I>j, for 0 < i < j — 1, is the normalizing transformation with 

respect to the periodic frequency w*. Here the periodic frequencies ujq, ... ,ujj are assumed to 
be known, Hj is already normalized with respect to uq, ... ,i^j-i and our aim is to explain 
the construction of the transformation ^j which will further normalize Hj with respect to 
ojj. The proof being technical, details will be given in AppendixO In ^3.31 a partial stability 
result (in the direction given by the linear span of coq, ... ,u!j, up to times of order 2™) will 
be easily deduced from the normal form Hamiltonian Hj o ^j. Using this, we will introduce 
a first version of the algorithm dichotomy in Proposition 13.21 either stability holds for an 
exponentially long interval of time and the algorithm stops, or a small drift in the action does 
appear in the orthocomplement of coq, ... ,LJj and the algorithm should move to the next step. 
We will also prove in Proposition 13.21 that if j = re — 1, then only the first alternative can be 
true. In 113.41 and 113.51 we will examine the situation where the second alternative holds true 
(so necessarily j < re — 2), and using the steepness property of h (in ^3.4p and then Dirichlet’s 
box principle (in ^3.51) we will prove how to pass from the step j to the step j + f. In 113.61 we 
summarize the work done in ^3.21 113.31 113.41 and 113.51 in Proposition 13.71 that gives one step 
of the algorithm, and which clearly shows that at some given step, either the algorithm stops 
or it yields the hypotheses that allow to apply it again. Finally, in TS.TI we conclude the proof 
of Theorem El which will follow easily from Proposition 13.71 by determining the parameters 
Q >1 and m > 1 in terms of our small parameter e. 

3.2 Normal form statement 

In this section we fix 0 < j < re —1, and we assume the existence of periodic vectors uq, ... ,ujj, 
with periods Tg,. .. Tj, which are linearly independent. For convenience we set W-i = 0 € M”. 
We define the complex and real vector space 

Aj := {u € C” I u • oj-i = V ■ ujQ = ■■■ = V ■ = 0}, Aj := Aj n M"", 

which are of complex (respectively real) dimension n — j. Then we consider three positive 
real numbers Vj, Sj and a point Zj G and we define the complex domain 

V3sj,35j(u) •= ^ I ^ l^(^) - Ikll < rj + 3^^} (3.1) 

where, for simplicity, the dependence on Vj is omitted. We will simply write |1 • || 3 sj-, 3 ^j. for 
the supremum norm for vector fields defined on V 3 s.^ 3 ^.(zj) and, for —1 < f < re — 1, we 
will denote by luji{z) := Ui ■ I{z) and its associated Hamiltonian vector held. With our 
convention, the function and its associated Hamiltonian vector held are identically zero. 

Given a real number 0 < e < 1 and an integer m > 1, we can dehne a set of Hamiltonians 
as follows. 
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Definition 3. The set NFj{uj-i ,..., coj-i, Zj, Sj,rj,^j, F, e, m), or for short NFj, consists of 
real-analytic Hamiltonians Hj defined on V 3 sj, 3 ^j{zj), and of the form 

Hj{z) := h{I{z)) + gj{z) + fj{z), z G V 3 sj, 3 ^j{zj) 

^h-.Vr^C, \\V^h{I)\\r<F, I{V3s,,3i,iZj))FVr, 

Il3sj,3?j < l|-^/jll3sj,3$j < j2^~^2~^e, 

^ ! 5j} ~ {^(^0 ’ 9j} ~ ■ ■ ■ ~ ) 5j} ~ 0- 

Now let us introduce another definition, taking into account the periodic frequency coj. 

Definition 4. The set NF j{uJo,... ,ujj, Zj, Sj,rj,f^j, F,e,m), or for short NFj, consists of 
real-analytic Hamiltonians Hj G NFj{ijj-i,... ,ujj-i,Zj,Sj,rj,f^j,F,e,m) which satisfy the 
following additional conditions: if we denote Ilj (respectively Hj) the orthogonal projection 
onto Kj (respectively Aj), then 

\\Iij'S/h{I{zj)) —ojjW = \\Hj'S/h{I{zj)) — ojjW < Sj (3-2) 


and 


Sj < {rj+ 2f^j)f^j, 2^21Q{rj + ?,f^j)mTj£<Sj, ^2{?,F^/n + l)f^N{rj-\-‘if^j)mTjSj<l. (3.3) 

The interest of the subset NFj C NFj is that if Hj G NFj, then np to a change of 
coordinates ^j (which is real-analytic, symplectic and close to identity), we get that HjO^j G 
NFj^i which will constitute a main ingredient in our algorithm. Here’s the precise statement. 

Proposition 3.1. Let Hj G NFj{uJo,... ,ujj, Zj, Sj,rj,f^j, F,£,m). Then there exist a real- 
analytic symplectic embedding 

: ^2sj,2^j{Zj) V3sj,3^j{Zj), <^j {V2sj,2^j{Zj)) D Vs^^^.{Zj), 

such that H^ := Hj o ^j = h + g^ + f(l~ with 

{L_i,5+} = = • • • = {4, 5 /} = 0, (3.4) 

and with the estimates 

\\Ng+\\2s„2i, < 2''+'e, ||^;+||2.„2€, < {j + 1)2^'2—e, (3.5) 

||^i-Id||2.,,2C, <2^^%£. (3.6) 

In particular, H^ G NFj^i{uj-i,... ,u]j, Zj+i, Sj+i,rj^i,f,j^i, F,£,m) given any choice of 
Sj+i, Vj+i, Cj-ei and Zj+i G for which the inclusion (z^+i) C V2sj,2ij{zj) 

holds true. 

The proof of Proposition 13.11 which is technical, is deferred to Appendix O The second 
part of Proposition 13.II follows easily from the first: if we define Hj+i := H(l~, gj^i := g^ and 
/j+i := then ()3.4p read 

{luj-ii 9j-ei} ~ = • • • = {li^j, gj-^-i} = 0 , 
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whereas the inclusion (z^+i) C V 2 sj, 2 ^j{zj) yields 

lO^Ssj+uSij+liZj+l)) ^ I{V2sj,2^jiZj)) c I{V3sj,3^j{Zj)) C Vr 

and, together with (|3.5I) . the estimates 

II^3j+iI|36,+i,3?j+i < \\Xgj^,\\2sj,2^j < 2^+^e, 

l|2^/,+ill3sj+i,3C,+i < \ \Xfj^^\\2sj,2ij < U + l)2^2"™e. 

This exactly means that -ffj+i = G NFj^i{uj-i ,..., Uj, Zjj^i, Sj+i, r^+i, F, e, m). 

3.3 Use of the normal form 

From now on, we will mainly work on the real domains 

V2sj,2ij (Zj) := V2sj,2(j (Zj) n 0 < j < n - 1. 

The normal form in Proposition 13.II is used to show that given a solution (t) of the Hamil¬ 
tonian system associated to = Hj o <hj, the curve I{z'^{t)) has a small variation in the 
direction spanned by wq, • • •, Wj, which is nothing but (the orthocomplement of Aj+i), 

for times |t| as large as the inverse of ||Ai^+|| 2 <jj, 2 ^j- It may well happen that in the direc¬ 
tion given by Aj+i, the curve I{z'^{t)) has also a small variation and hence I{zj{t)) where 
Zj{t) = <hj(z^(t)), has small variation, which yields our confinement result. But if not, that 
is if there is a faster deviation of I{z^{t)) from l{z^{0)), this has to occur in the direction 
given by Aj+i. Here is a precise statement. 

Proposition 3.2. For 0 < j < n — 1, let Hj G NF j{u}o,... ,ujj, Zj, Sj,rj,^j, F, e,m) and 
: V 2 sj, 2 ^jizj) —)• Vssj, 3 ^j{zj) given by Provosition \8.1\. and let z'^{t) be the forward solution 
of the Hamiltonian H^ = Hj o ^j starting at z'^ := ^J^{zj). If 0 < j < n — 2, and if we 
define 

ij := (r, + + l)-i2-^s-i2-, (3.7) 

then we have the following dichotomy: 

(1) either z^(t) G Vs.^^.{zj) for 0 < t < ij, 

(2) or there exists a positive time tl < ij such that 

-Hh)\ = ■Si/4 

and, for 0 < t < t^, 

z^{t) G Vs^,^.{zj), \I{z+{t))-I{zj)\ < Sj/4, \Uf^^{I{z+{t))-I{zj))\ < sj^e. (3.8) 

If j = IT'~ if we define 

in-i := {rn-i + Sn-ie~^2‘^, (3.9) 

then z:^_^{t) G (^^n-i), forO <t< in-i- 
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Proof. First observe that since the image of contains Vs.^^.{zj), it contains Zj and so 
z'^ = ^J^{zj) is well-defined, and we have, using (I3.6h and the hrst two inequalities of (I3..SD . 

\\z+ -Zj\\ = ||z+ - ^j{z+)\\ < < s,(108(r,- + 3?,))-^ < + 2Cj)(108(r,- + 3^,))-^ 

which easily implies 


114 - ZjW < ?i/108, \I{z+) - I{zj)\ < Sj/108. (3.10) 

Observe also that since zj is real and Hj and are reals, the forward solution z^ (t) is real. 

We first consider the case 0 < j < n — 2. Using (13.101) . we can now define G (0, -|-oo] to 
be the time of first exit of z^ (t) from i^j)- We claim that the dichotomy of the statement 
is implied by the following trivial dichotomy: either tj < tj or tj < tj. 

Indeed, in the first case, one obviously have 


z+ (t) G W. (zj), 0<t<tj. 

In the second case, either \I{z'^(tj)) — I{zj)\ = Sj or ||z)''(fp|| = rj+^j. But since the solution 
is real, the second possibility implies that 


|/( 3 +(t'))l = l/2||2+(tj)||2 = l/ 2 (rj + 4 j) + 

while 

\I{zj)\ = l/2\\z,\\^<l/2r] 
and therefore, using the first inequality of (13.3|) . we obtain 


|/{z+(t')) - I(z,)\ > \I{z*[f‘))\ - |/{z,)| > {,(l/2f, + r,) > s,V4. 

So, whether \I{z~j{t^)) — I{zj)\ = Sj or l|2:j'(tj)|| = Vj -\- there exists a positive time 
— h such that 

-Hzj)\ = Sj/4: 

and 

\I{z+{t))-I{zj)\<sj/4, 0<t<t+. 

Since <tj, G for 0 < f < It remains to show that 

\Uf^^{I{z^{t)) - I{Zj))\ < sj^£, 0<t<t+. 

Since h is integrable, for 0 < s < f < 


fL 

dt 


I{zj{s)) 


{/, H+}{z+{s)) := {{h,H+}{z+{s )),..., {4, H+}{z+{s))) 
{/, h + g+ + f^}{z+{s)) = {/,5+ + //}(<(s)). 


Then, for any z G Vs.^^-{zj) and any — 1 < f < j, using (|3.4p we obtain 


=0Ji • {I,5(+}(z) = 0 
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which implies that {I,g^}{z) G Aj+i. Therefore 

ni+i = {^’//}(<(^)) 

and hence 

and, using the second inequality of (|3.5h and the fact that 

< ii^;wiiiix/(2;w)ii 

we obtain 

- I{zj))\ < t{rj + Cj)\\Xf+\\ 2 sj, 2 ^j < t{rj + ^j){j + 1)2^2-^£. 

Since t'j' < tj < ij = {rj + + 1)~^2~^sj^2^, we thus obtain 

\Uf^^{I{z^{t)) - I{zj))\ < sj^e, 0 < f < f+, 

which concludes the proof for the case 0 < j < n — 2. 

Now for the case j = n — 1, we have A„ = {0} and so = M"', hence for t < tn-i = 
(vn-i + repeating the last argument we get 

|^(^n_l(i)) - I{Zn-l)\ < Sn-l, 0 <t < tn-l- 

As z^_i{t) is real, this implies, using also the first inequality of (I.S.3p . that for 0 < f < fn-ij 
Il4_,(f)l|2 = 2|/(4_,(f))| 

< 2|/(4_i(f))-/(z„_i)| + 2|/(z„_i)| 

^ ‘2Sn-l + \\Zn-l{t)\'^ 

^ 2(r^_l + 2^n—l)Cn—l + ^n—1 

< (r„_i + 2en-i)^ 

so2:^_^(f) G V's^_i,25„_i(-2n-i) for 0 <t< tn-i, and this concludes the proof of the proposition. 

□ 


3.4 Use of the steepness property 

Let us start by giving a geometric interpretation of the steepness property, as its definition is 
quite abstract. Assume that h is steep on some domain D, and consider a curve 7 : [0,1] —?■ M” 
which takes values in Xf] D, where A is a proper affine subspace of M”. It may happen that 
V/ia( 7(0)) = 0 (this is the case if 7(0) is a resonant point for h, that is, if k ■ V/i(7(0)) for 
some non-zero integer vector k G 1A: then V/ia( 7(0)) = 0 where A is the real space generated 
by such integer vectors k). If this happens, the steepness property ensures that, for some 
time 0 < f < 1, X/hx{'y(t)) 7^ 0 (informally, in terms of resonances, this means that we do not 
have “accumulation of resonances”). Moreover, the longer is the length of the curve 7, the 
farther away from zero is the vector Here’s a quantitative statement, which is due 

to Nekhoroshev. 
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Lemma 3.3 (Nekhoroshev). Let h he a function whieh is (r, At, C, <5, and 

such that 

\\V'^h{I)\\r < F. 

Let 7 : [0,t~^] —)• M” be a continuous curve, A an affine subspace o/M” of dimension I, where 
1 < I < n — 1, and d a positive real number. Assume that 

(i) for all t S 'y{t) S A; 

(u) for all t G [0,t+], ||7(0) - 7 (t)|| < d and ||7(0) - 7(t+)|| = d; 

(Hi) the ball {L G M"" | — 7(0)|| < d} is contained in Dj.; 

(iv) d < min{(5, (3 F)“^k, 2(5/t(4C')“^)^/^’'}, 
then there exists a time t G [0, t~^] such that 

||nAV/i(7(t))|| >c/5(d/2r, 

where A is the vector space associated to and IIa the orthogonal projection onto A. 

This is a special case of the lemma on “almost plane curves” of Nekhoroshev, stated in 
|Nek77| and proved in [Nek79] (our case corresponds to “plane curves”). 

Now assume that Alternative (2) of Proposition 13.21 holds true, and let 7j(t) := + 

Ilj^i{I{z^{t)) — L{z^)) fort G [0,t^]. Since this curve takes values in a proper affine subspace, 
the following proposition is a simple consequence of Proposition 13.21 and Lemma 13.31 


Proposition 3.4. For 0 < j < n — 1, let Hj G NFj{u}o,... ,ujj, Zj, Sj,rj,f,j, F) and : 
V 2 sj, 2 ^jizj) —>■ V 3 sj, 3 ^j{zj) given bu Pronosition HOI and let z^ ft) be the forward solution of the 
Hamiltonian H'j' = HjO^j starting at z'^ = ^~^{zj). Assume that h is (r. At, C, 6, {pi)i=i,...^n-i)- 
steep. Then we have the following dichotomy for j < n — 2: 

(1) either z^ (t) G (zj) for0<t< ij, 

(2) or there exists a time ij <ti'< ij such that, setting Jjiij) '■= I{zli) + Ilj^i{I{zli{ij)) — 


I{zj)), then 


provided that 


|nj+iV/i(7j(tj))|| > fijs 


Pn—j—l 




fe < spS, 

\sj < 8min{(5, {3F)-^k,2{5k{4.C)-^)^/p^-^-^}. 
Iff = n-l, then 2+_i(t) G for0<t< t^-i. 


(3.11) 

(3.12) 


Proof. We only have to consider the case j < n — 2, and we have to prove that Alternative 
(2) of Proposition 13.21 implies Alternative (2) of the above proposition. So we assume the 
existence of a positive time tf < i.j such that 


1^(4(4)) -P4\ = ■Si/4 
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and, for 0 < i < 

^ - -^(^i)l < «i/4, |nj-+iU(^/(i)) - ^(^i))l < sj^e- 

Hence, using the first inequality of (|3.12l) . 

|ni+i(/(^/(i+))-/(^i))| > \I{z-^ {tp) - I{zj)\-\Yif^^{I{z-^ {t)) - I{zj)\ > Sj/A-sJ^e > Sj/8 
and in particular 

\\Uj+,{I{z+{t+))-I{zjm>s,/8. 

Therefore we can certainly find a positive time ij < t'j' such that 

\\U,+i{I{z+ii+))-Iizjm=Sj/8 

and 

\\Uj+i{I{z+it)) - I{zjm < Sj/8, 0<t< i+. (3.13) 

Now we want to apply Lemma [H31 to the curve 7 j(t) = ) — I{z^)), for 

t G [0,t^], with d := Sj/8 and with the affine subspace Aj+i := I{z^) + Aj+i which has 
dimension n — j — 1. The assumptions (i) and (ii) of Lemma (13.311 are trivially satisfied, as 

ijitp - 7i(o) = nj+i(/(z+(t+) - i{zj)). 

Then (in) holds true since, by definition of Hj, we have I{V 2 sj, 2 ^jizj)) C 3 ^^. (z^)) C Dr- 

Eventually, the second inequality of (I3.12|) clearly implies {iv) therefore Lemma 13.31 can be 
applied, and there exists a time ij G [0,t)^] such that 

||n,+iV/i(7i(t,))ll > 5-^C{d/2r-^-^ = 

This concludes the proof. □ 

3.5 Use of periodic approximations 

Let us Hrst state the following simple consequence of Dirichlet’s theorem on approximation 
of real vectors by rational vectors. 

Lemma 3.5. Let v G M” \ {0}, and Q > 1 a real number. Then there exists a T-periodic 
vector uj G M” \ {0} such that 

||u — w|| < Vn— 1{TQ)~^, Ibir^ <T< y/n\\v\\~^. 

Proof. Fix <5 > 1. Up to a re-ordering of its component, we can write v = |?;|(± 1 , 3 :) for 
some X G and by Dirichlet’s approximation theorem, there exists a rational vector 

p/(? G such that 

\qx-p\<Q~^, l<q<Q'^~^. 

The vector u) = |u|(±l,p/g) G M” is then T-periodic, for T = \v\~^q, and we have 
llr-wll <T-^\\qx-p\\, \v\-^<T< \v\-^Q'^-^ 

which implies 

||u — cull < y/n — ll^ll”^ < T < y/n\\v\\~'^Q'^~^ 

and this was the statement to prove. □ 
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Now assume that the conclusion of Alternative 2 of Proposition 13.41 holds true, so the 

vector nj_|_iV/i( 7 j(tj)) is non-zero, where 7j(tj) = — I{z^)). By 

Lemma [3.51 this non-zero vector can be approximated by a periodic vector Wj+i, and it will be 
easy to ensure that this new periodic vector is linearly independent from wq, ■ ■ ■ ,tdj (as Wj+i 
is close to Ilj+i'Vh{'jj{tj)), the latter being, obviously, linearly independent from ujq, ... ,ujj 
as it is orthogonal to them). Moreover, as ^j{ij) is close to I{z^{ij)), setting zj^i := z^{ij), 
the vector (Uj+i is also an approximation of nj+iVh(/( 2 :j+i)). This leads to the following 
proposition. 

Proposition 3.6. For 0 < j < n — 1, let Hj S NF j{u}o,... ,LJj, Zj, Sj,rj,^j, F, e,m) and 

given by Proposition and let z^{t) be the forward so¬ 

lution of the Hamiltonian Hf = Hj o <!>,• starting at zf = ^~^(zj). Assume that h is 
(r, K, C, 6, {pi)i=i,...,n-i)-steep. Then, for j < n — 2, we have the following dichotomy: 

(1) either z^(t) G Vs.^^.{zj) for 0 < t < ij, 

(2) or, given some positive parameter Q > 1, there exists a Tj+i -periodic vector Wj+i G 
M" \ {0}, linearly independent from uiq, ... ,u:j, with the estimate 

{F'y/flsj)-^ < Tj+i < F' := max{l, F], (3.14) 

and a time tj < < ij such that if we define 

Zj-j-i '.= Zj {ij), ■Sj-i-i := 2-\/n — l{Tj^iQ) , Vj^i := rj -|- f^j, := 'Cj/3 

then it holds that 

^3sj+i,3^j+l{Zj+l) QV2sj,2^j{Zj) (3.15) 

\\Uj+iVh{I{zj+i)) - uij+iW < Sj+i, (3.16) 

z^it)eVs,,^^izj), 0<t<ij, (3.17) 

provided that 

e < sj/S, 

^ Sj < 8min {J, (3 F)- ^k,2(5k(4C)-^)Vp-^-i}, ^ 

Q > SF'■sjn{n — 1), 

e < {2^/nF)~'^SjSj+i. 

If 3 =n-l, then z:^_^{t) G 14„_i,€„_i(2:n-i), for0<t< in-i- 

Proof. Since (13.181) implies in particular (13.121) . it is enough to prove that Alternative (2) of 
Proposition 13.41 implies Alternative (2) of the above proposition. So we assume the existence 
of a time ij < < ij such that 

||n,+iV/i(7,(t,))|| > = 5-'C16-^’—1, 

where 'Jjiij) = I{z^) + Hj-eiili^^(ij)) ~ H^j'))- Let us define Vj+i := Ilj^i'Vh{'-fj{ij)). We 
have 

W^jitj) -I{zj)\\ < \\Uj+i{I{z+{ij)) - /(z+))|| -h 11/(2+) -I{zj)\\ < Sj/8 + y/nsj/108, 
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where we used the estimate (j3.13l) (as ij < and the estimate (I3.10p . 
implies 


\hj{tj) - Iizj)\\ < {^/n-l)sj 


Since n > 2 this 


which implies 


Ibi+i - nj+iVh(/(zj))|| < F{^/n-l)sj. 

Now recall that by definition of zj, we have 

\\IljVh{I(zj)) — iOjW < Sj 

and since Uj G and A^+i C Aj, Ilj+iWj = 0 and Hj+i = Ilj+iIIj, and therefore 

||n,-+iVh(/(z,))ll = \\U,+,{U,{Vh{I{zj)))-u;j)\\ < s, 


which implies that 


||n,-+i|| < ||n,+i -n,+iVh(/(z,))|| + ||n,+iVh(/(2,-))ll < F'V^SJ. 


We just proved that 

< ||nj+i|| < F'y/nsj. (3.19) 

Now, for <5 > 1, we apply Lemma [3.51 to Uj+i: there exists a Tj+i-periodic vector Wj+i G 
M" \ {0} such that 

\\vj+i - ujj+iW <Vn- l{Tj+iQ)-^, lluj+ill"^ <^i+i < Vn\\vj+i\\~^Q'^~^. 

Using (13.191) . this implies 

{F'y^Sj)-^ < \\vj+i\\-^ < Tj+i < y^\\vj+i\\-^Q^-^ < (3.20) 

and also, nsing the lower bound on Fj+i, 

\\vj+i -ujj+i\\ < y/n- l{Tj+iQ)~^ < y/n - l\\vj+i\\Q~^ < F'y/n - ly/nsjQ~^. (3.21) 

Let ns prove that Wj+i is linearly independent from cuq, ... ,ujj, that is ujj+i does not belong 
to To do this, it is enough to prove that if v is an arbitrary vector in then 

\^jJj+i ■ v\ < ||t<;j_(_i||||n||: indeed, otherwise, letting v = Wj+i, one wonld get a contradiction. 
On the one hand, we have 

\LOj+i ■ v\ = \{ojj+i - Uj+i) • n| < \\vj+i - Wj+illllnll < y/n - lQ“^||nj+i||||n|| 
where we used the fact G Aj+i and (I3.2ip . while, on the other hand, 

||wj+i||||n|| > (||nj+i||-||nj+i-a;j+i||)||n|| > {1-y/n - lQ-^)\\vj+i\\\\v\\ > y/n - lQ-^\\vj+i\\\\v 

where we nsed the third inequality of (I3.18p and (|3.2ip . These last two inequalities imply that 
Icuj+i • v\ < ||a;j+i|l||n|| for an arbitrary vector v G Aj:)_^, and so Wj+i is linearly independent 
from ujQ ,... , ojj. 

Next we dehne 

Zj+i := z^iij), 
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and observe that, since tj < by (|3.8I) . Zj+i G but also G Vsj^^^(zj) for 

0 < t < ij, which justifies ()3.17p . Moreover, still from (|3.8h . 

\I{zj+i) - 7j(tj)| = |n^+i(/(zj+i) -I{zj))\ < sjie 


so 

WHzj+i) - 7j(ij)ll < VnsJ^e 

and hence 

\\Uj+iVh{I{zj+i)) -Vj+i\\ < Fy/ns~^e. 

Therefore, using the first inequality of ()3.21h . the definition of Sj+i and the last inequality 
of (IM) . 

\\Uj+iVh{I{zj+i)) - Uj+iW < \\Uj+iVh{I{zj+i)) -Vj+i\\ + \\vj+i - ujj+i\\ 

< F^/ns'J^e + Sj+i/2 < Sj+i, 

which proves (j3.16h . It remains to check (I3.15p . so let us hx z G First, we 

have 


I{z) - I{zj) = I{z) - I{zj+i) + I(zj+i) - I(zj) =I{z)- lizj+i) + I{z^{ij)) - I{zj) G Kj 

since I{z) — I{zjj^i) G Aj+i C Kj and I{z^{ij)) — I{zj) G Aj. Then, 

\I{z) - I{zj)\ < \I{z) - I{zj+i)\ + \I{zj+i) - I{zj)\ < 3sj+i + Sj/4 < 2sj 

provided that Sj+i < 7sj/12: but this inequality (in fact, the stronger inequality Sj+i < Sj/2) 
follows from the definition of Sj+i, the third inequality of (I3.18p and (|3.21l) . Eventually, 

11^11 ^ A ~ A A ~ A 

and so we showed that (zj+i) C V 2 sj, 2 ^jizj), which concludes the proof. □ 

3.6 One step of the algorithm 

As a straightforward application of Proposition 13.11 Proposition 13.21 and Proposition 13.61 we 
now describe formally one step of the algorithm that will eventually lead to the proof of 
Theorem [Dl 

Proposition 3.7. For 0 < j < n — 1, let Hj G NFj{uJo,... ,ujj, Zj, Sj,rj,^j, F) and : 
V 2 sj, 2 ^j{^j) '^ 3 sj, 3 ^jizj) given by Provosition \d.li and let z^{t) be the forward solution of the 

Hamiltonian = Hjo^j starting at z^ = ^J^{zj). Assume that h is (r, n, C, 6, (pO/=i,...,n-i)- 
steep. Then, for 0 < j < n — 2, we have the following dichotomy: 

(1) either z^{t) G Vs.^^.{zj) for 0 < t < ij, 

(2) or, given a real number Q > 1, there exists a Tj^i-periodic vector Wj+i G M” \ {0}, 
linearly independent from ujq, ... ,ujj, with the estimate 

< T,+i < F' = max{l, F}, (3.22) 


30 
















and there exists a time ij < tj such that, 


Zj-\-i — {ij) G 'Sj+i — 2\/re l{Tjj^iQ) , ^j+i — 'fj + ‘^j+i — 

such that H'i G NFjj^i{ijjQ, ..., cjj+i, ^j+i, Sj+i, r^+i, F, e, m) and 

z-^{t)eVs,,^,{zj), Q<t<ij, (3.23) 

provided that 

'Sj < 8min{(5, (3^)-^^, 2(5 k(4(:7)-1)1/P"->-i}, 

Q > SF'■sjn{n — 1), 

< e < {2y/nF')~^SjSjj,.i, (3.24) 

2^+^2lQ{rj+i + 2,ij+i)mTjj^ie < Sj+i, 

72{‘iFy/n + l)^“j!i(rj+i + ?,ij+i)mTj+iSj+i < 1. 

If j = n-l, then 2;[_i(t) G (^n-i) forO <t< in-i- 

Proof. The case j = re — 1 follows directly from the case j = re — 1 of Proposition [321 Then, we 
claim that the inequalities (13.241) imply the inequalities (13.181) and the inequalities (13.3p (with 
j replaced by j’ + 1). Assuming this claim, and using Proposition 13.61 we have the inclusion 
of the complex domains of (j3.15l) . and therefore using the second part of the statement of 
Proposition 13. 11 we can assert that G NFj^i. Moreover, in view of (j3.161) . and since (|3.3I) 

is satisfied (with j replaced by j + 1), we eventually obtain that G NFj^i, while (j3.23l) 
is nothing but (13.171) . 

It remains to prove the claim. To do this, observe that (I3.24P obviously implies (13.181) 
and (j3.3l) . except for the following two inequalities: 

e < s^/8, Sj+i < (rj+i + 2^j+i)^j+i. (3.25) 

But using the third inequality of (13.241) and the fact that Flj G NFj , we know that 

e < {2^/nF')~^SjSj+l, Sj < (rj + 2^j)^j. (3.26) 

Then, using the second inequality of (13.241) . one easily check that Sj+i < Sj/A, and this, 
together with (I3.26p . imply (I3.25p . and the proof is over. □ 

3.7 Proof of Nekhoroshev exponential stability 

We can finally give the proof of Theorem [Dj Recall that we are given a Hamiltonian H as 
in d*]), which is dehned on Br, and of the form 

H{z) = h{I{z)) + f{z), h:Vr^C, f : Br ^ C 

and that pi. 81) holds true, that is 

\\Vh\\r<E, \\V^h\\r<F, \\Xf\\r<e. 
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We already defined F' = max{l,F}. Recall also that h is (r,n,C,6, 

us now define additional parameters; for any 0 < j < n — 1 and 0 < k < j, we set 

■= n “i •= 

n—j<i<n—j-\-k—l 0<i<k 

with the convention that the product over the empty set is one, that is, ir^ = 1. Observe in 
particular that 

Oq = 1, a\ = l+pn-i, 

and at the other extreme, 

a”l2 = 1 + P2 + P2P3 H-h P2P3 ■ ■ ■ Pn-1 = a', 

a^z\ = l+Pl+ PlP 2 H-h PlP2 ■ ■ ■ Pn-l = a. 

For 0 < j < n — 2, recalling that the numbers pj have been dehned in (|3.11H . we dehne 

V ■= min {(3 F)-^k,2(5k(4C)"^)^/P"-^-i}, i 2 j := 

0<j<n—2 

2 = 0 

The proof of Theorem iDl will be a consequence of the following proposition. 

Proposition 3.8. Let H{z) = h{I{z)) + f{z) be as in d*]) satisfying (|1.8I) . such that h 
is {r,K,C,6,{pi)i=i^,,,^n-i)-steep. Let zq be an arbitrary point in Bj ./2 and z{t) the forward 
solution of H starting at zq. Given an integer m > 1 and a real number Q > 1, we have 

\L{z{t)) — I{zo)\ < s := 3Ey/n — 1Q~^, 0 < t < t := 3{2rE\/n — 1)~^Q2'^, 


provided that: 


'Q > (5r2)-l36£;^/n^, 

Q > E\/n — 1(8?/)“^, 

Q > .FVn^(8(5)-\ 

Q > SF' n(n — 1), 

2n-i27(3 _j_ 3-"-+i)7'(rj — L)~°-n°'y/n — e < 1, 

> m216(3” + l)(3Fv^+ l)v^nr^. 


(C) 


Let us first prove this proposition. The fact that this proposition implies Theorem [D] 
simply follows from a suitable choice of m and Q (in terms of our given parameters) and will 
be detailed later. 


Proof of Proposition I,'?. <91 The proof follows from an algorithm whose inductive step is given 
by Proposition 13.71 But hrst we need to initiate the algorithm. By assumptions we have 

K<\\^h{I{zom<E 

and so we can apply Lemma 1331 to vq := Vh{L{zo)): there exists a Tg-periodic vector ujq G 
M" \ {0} such that 

ll^^o - c^oll < V^i^{ToQ)-\ E-^ < To < (3.27) 
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We define 


Ho:=H, So ■■= Vn^{ToQ) \ ro := r/2, := J’o/S = r/6, 

and observe that Hq G A^Fo(a;_i, zq, sq, ro, F, £; "i). Indeed, Aq = C”, ro + 3,^o = r so that 
V 35 o, 3 Co(^o) C Br, and we can write Hq = h + f = h + go + fo, with 50 := / and /o := 0, 

ll^5oll3ro,35o — ll^/l|s — *^5 

as the requirement {lui_i,go} = 0 is void since U-i = 0. In fact, using the first inequality 
of (|3.27p and assuming that 

'Q > (5r2)-i36£;Vn^, 

< 216ny/n — 1 ^K~'^rmQ‘^^~^e < 1, (CO) 

Q > m216.3{3Fy/n + l)\/?r — 1, 

one easily check that, using the dehnitions of sq, rg, ^0 (which gives in particular rg + 3^o = r 
and (ro+3^o)'^o~^ = 6) and the second estimate of (I3.27p . that Hq G NFo{ljJo, zq, soj ro, ^ 0 , F, e, m). 
So Proposition 13.71 can be applied. 

If Alternative (1) of Proposition l3.7l holds true, the solution z^ {t) of Hq = Hqo^q satisfies 
z^{t) G I4o,?o(^o) for 0 < t < to- As d>o sends l 2 so, 2 ?o(^o) into P3so,3Co(2o) and t < to, then 
<I>o(2(('(t)) = zo{t) = z{t) satisfies in particular 

\I{z{t)) — I{zo)\ < 3so < s, 0 < t < t, (3.28) 

the proposition is proved and the algorithm stops. 

If Alternative (2) of Proposition 13.71 holds true, then there exist a Ti-periodic vector 
coi G M”' \ {0}, linearly independent from ujq with the estimate 

(F'V^so)-^ <Ti< (3-29) 

and 

Zi = z^(io) € Br^, Si = 2y/n- l{TiQ)~^, ri=ro + ^o, 6 = ^o/3 

such that G A'Fi(ti;o, wi, ^i, si, ri,^ 1 , F,e, m) and 

Zoit) &Vso,^o{zo), 0<t<io, (3.30) 


provided that 


Q > E\/n — 1(8??) 

Q > Ey/^r^{8S)-\ 

Q > 8F'^Jn{n — 1), 

180r(n — l)““in“i-v/n — < 1, 

Q > m216.10(3FY/n + l)y/n — 1. 


(Cl) 


Indeed, using the definitions of so, si, rg, ri, ^g, (in particular, we use the facts that 
Si < sg. Si > y/n — l{TiQ)~^, ri + 3^1 = 5r/6 and (ri + 3^i)^j"^ = 15) and the estimate (j3.27h 
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and (13.291) on respectively Tq and Ti, one can check that (ICll) imply (13.241) for j = 1. Setting 
H, := H+ e NFi{iOo,LOi, zi, F,e,m), we can apply Proposition 13.71 again. 

If Alternative (1) holds true, then the solution zf{t) of = Fli o = Hq o 4 >i = 
Hoo ^qo $1 starting at zf = satishes zf{t) G for 0 < t < ii. As <I>i sends 

V2si,2a(^i) ii^to P35 i,36(2:i)> then 4 >i(zj''(t)) = zi{t) belongs to 143^1,3^1(2:1) for 0 < t < h- 
By (|3 15p . 143si,3^i(2;i) is contained in 142so,2?o(^o), and as t < ii, zi{t) belongs to V2so,2^o{zo) 
for 0 < t < t. Now observe that since 2:1 = 2()'(to), by uniqueness of the solutions associated 
to the system dehned by Fli = Fl^, we have the equality 2:1 (t) = 2;Q"(t + to) as long as the 
solution is defined. Using this equality, what we have proved is that 

z^(t) G 142 ^ 0 , 2 ^ 0 (^ 0 ), k <t <io + i 
But recall that from (j3.3np . we know that 

Zo{t)eVso,^oi^o), 0 < t < to, 
and therefore, since t < to +1, we have in particular 

z^ (t) G 142^0,2So(^o), 0 < t < t. 

As before, using this and the fact that $0 sends 142so,2^o('^o) into 143 so,3?o('^o) we also arrive at 
the estimate (I3.28|) . 

If Alternative (2) holds true, then the algorithm continues. To apply Proposition 13.71 at a 
step j, for 1 < j < n — 1, it is sufficient to check that 


Q > Ey/n — 1(87?) 

Q > p;^^r4ri(85)-i, 

Q > 8F'y/n{n — 1), 

O /— TP> /— 1 /-i —— —1 —1 —(7rl+7rl~l)^n(al+al~l) ^ -t 

2y/nl< y/n ^ ^-^y/n — 1 ^ ^ V- ^ ^ Q ^ < 1, 


(Cj) 


j-i 

2^27(3 + 3-J>(n - 1)““^ 

Q > m216(3J+i + l)(3FVn + l)Vn44I. 


-\2na^- — 1 


e< 1, 


Indeed, ( |CjP implies (13.241) . using the definitions of Sj, rj and for 0 < i < j (which imply 
in particular that the Si are decreasing, Si > y/n — l{TiQ)~^, ri + 3.^i = r(3 + 3“*)/4 and 
{ri + 3.^1)^“^ = 3(3*41 _|_ i)/2), and the estimates on the period Tj that one obtains at each 
step using (|3.22l) . To conclude, just observe that the conditions ([U]) imply the conditions (ICOIl 
and (Cj) for any 1 < j < n — 1. For j = n — 1, there is only one possibility in Proposition (13. 7p . 
the algorithm stops and the statement is proved. This ends the proof. □ 

Proof of Theorem\^ We just need to choose m and Q in Proposition 
given parameters. First we choose m in terms of Q as follows: 


in terms of our 


where 


m := [biQ], bi = (216(3Fv*n + 1)(3"' + l)y/n — 1) ^ 
denotes the integer part. Using this choice, the conditions (jC]) are implied by 


Q>62, Q>M"\ Q>bAr-\ r65Q^^“e<l, beQ^^'^+‘^'h < 1, (3.31) 
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where 


62 := max{8F'Y^n(n — l),Ey/n — 1(877)“^, 6^^} 
63 := 

64 := 5“^36£\/n — 1 
65 := 2”-i27(3 + 3-"+i)n“(n - 

fte := 2-vAiF'-v/n“^“ Vn - 1 
Then we choose Q follows: 

Q := {b^rey^ 

and observe that () 3 . 3 ip becomes 


2a a + a 
a —a' T Cl- 

^6 °5 


re < re < re < , re < 6g br’^'r—' 

With these choices of m and Q, since m > biQ — 1 we have 

1 I 

s = 3 Ey/n — 16 |"“ (re)^ 

and 


t > 3 {AE\/n — 1 ) ^(re) 2na exp ^(In 


2)6165 ^”‘°-r~^{re)~'^ 


so if we dehne 


ci:=65-V”“, C2 := 65-I63-2-, c ,-.= := b , 


2a a + a' 

a — a^ a —a' 


6 h 


and 


C5 


:= lE+n — 165"-“ , C6 := 3 ( 4 £i\/re — 1 ) ^65 cy := (1112)6165 


1 

2na 


we eventually obtain that if 


re < min |ci, 026^”“, C3r'^”“, C4r“-“' | 


then 


( 3 . 32 ) 


( 3 . 33 ) 


|/(z(t)) —/(^;o)| < C5(re)2na j 0 < t < cgr ^(re) 2na exp ^cyr ^(re) 2na^ . 

This proves the statement for positive times, but for negative times, the proof is of course the 
same, so this concludes the proof. □ 
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A Proof of generic steepness 


The aim of this section is to give the proof of Theorem 12.11 The latter will be an immediate 
consequence of ProDositions lA.2l and lA..‘il below. We shall use in the proof of these propositions 
basic results concerning semi-algebraic subsets; for proofs and more information we refer to 
|BCR98| . Our main ingredient to prove Theorem I2.1I is a result of Nekhoroshev on stably 
expanding polynomials that we will now state. 

Let us first recall that P{n, m) denotes the space of polynomials of degree m in n variables 
with real coefficients, and P 2 {n,m) the subspace of P{n,m) consisting of polynomials with 
vanishing homogeneous parts of order zero and one. The following definition, which is related 
to the definition of stably steep polynomials, is due to Nekhoroshev f [Nek73| L 

Definition 5. Let 1 < I < n — 1. A polynomial Qq G P 2 {l,m) is called stably expanding if 
there exist a neighborhood Ui ofQo in P 2 {l,m) and positive constants C[,5[ such that for any 
Q G Ui, the inequality 

max min ||VQ(y)|| > C/r"' 
o<u<i ||y||=f? 

holds true for all 0 < < 6[. 

The set of stably expanding polynomials in P 2 {l,m) will be denoted by SE{l,m). 

Theorem A.l (Nekhoroshev). Let 1 < I < n — 1. The complement of SE{l,m) in P 2 {l,m) 
is contained in a closed semi-algebraic subset T,{l,m) of codimension [m/2]. 

Let us denote by L{n, 1) the space of rectangular matrices with n rows and I columns, with 
real coefficients, and by Li(n, 1) the open subset of L[n, 1) consisting of matrices of maximal 
rank. Any A G L{n,l) induces a linear map A : —)• M”, hence given P G P{n,m), we can 
define Pa G P{l,m) by setting Pa{x) = P{Ax), x G Moreover, if P G P 2 {n,m), then 
Pa G P 2 {l,rn). Let us define the set 

0(/,n,mo) = {{P,A,Q) G P 2 {n,mo) x Li{n,l) x E{l,mo) \ Pa = Q}. 

Then we define T{l,n,mo) to be the projection of 0(Z,n, mo) on the first factor P 2 {n,mQ), 
and finally 

n—1 

T(n,mo) = IJ T(Z,n,mo). 

i=i 

Theorem 12.11 is a straightforward consequence of the following two properties of the set 
T(n,mo). 

Proposition A.2 . The set T(n,mo) is a semi-algebraic subset of P 2 {n,mo) of codimension 
at least one. 

Proposition A. 3. The complement of SS{n,mo) in P 2 {n,mo) is contained in T(n,mo). 

The second proposition is true for any m > 2 and not just for m = mo, but this will not 
be needed. 

Let us now give the proof of Proposition lA.21 and Proposition ! A. 31 following the arguments 
in |Nek73| . 
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Proof of Provosition \A.^ The set P 2 {n, mo) is a real vector space hence it is algebraic, Li{n, 1) 
is obviously an algebraic subset of L{n,l) whereas, by Theorem lA.ll is a semi- 

algebraic subset of P 2 {l,mo). Moreover, for {P,A,Q) G P 2 {n,mo) x Li{n,l) x the 

equality Pa = Q corresponds to a system of algebraic equations in the coefficients of P, A and 
Q. This implies that Q{l,n,mo) is a semi-algebraic subset of P 2 {n,mo) x L{n,l) x P 2 {l,mo). 
Now since the projection of a semi-algebraic subset is a semi-algebraic subset, T{l,n,mo) 
is a semi-algebraic subset of P 2 {n,mQ). Then, as a finite union of semi-algebraic subsets is 
semi-algebraic, T(n,mo) is a semi-algebraic subset of P 2 {n,mo). We need to prove that the 
codimension of T(n,mo) in P 2 {n,mQ) is at least one; to do this it is sufficient to prove that 
the codimension of T(^,n,mo) in P 2 {n,mo) is at least one for any 1 < ^ < re — 1. So let us 
fix 1 < / < re — 1. Given {A,Q) G L{n,l) x P 2 {l,mo), we define 0A,Q(^,n,mo) to be the 
intersection of Q{l,n,mo) with the set 

{{P',A',Q') G P 2 {n,mo) x L{n,l) x P 2 {l,mo) \ A' = A, Q' = Q}. 

If {A,Q) G Li(re,/) X S(/,mo), it is easy to see that 

dim0^^Q(/, re, mo) = dimP 2 {n,mo) — dimP 2 {l,mo) 

and therefore 

dim0(Z,re, mo) = dim0yi^Q(/, re, mo) + dimLi(re, ^)-|- dimS(^, mo) 

= dimP 2 {n,mo) — dimP 2 {l,mo) + dimLi{n,l) -|- dimS(?, mo) 

= dimP 2 (^) iT^o) + dimLi(re, 1) — codimS(/, mo) 

= dimP 2 (?T') n^o) +nl — [mo/2] 

where in the last equality we used the fact that dimLi(re,/) = dimL(re,^) = nl and Theo- 
rem lA.ll Now given P G P 2 {n, rno), we define 0p(/, re, mo) to be the intersection of 0(/, re, mo) 
with the set 

{{P',A',Q') e P 2 in,mo) X L{n,l) X P 2 {l,mo) \ P' = P], 

Recall that if GL[l) denotes the group of square invertible matrix of size I, with real coeffi¬ 
cients, then GL{1) acts freely on Ti(re, 1) (the quotient space is nothing but the Grassmannian 
G{l,n), that is, the space of all /-dimensional subspaces of M"). It is then easy to see that 
GL{1) acts freely on Qp{l, re, mo), therefore the dimension of an orbit of this action equals the 
dimension of GL{1), which is P, and hence, 

dim0p(/, re, mo) = P. 

Since T(/,re,mo) is the projection of Q{l,n,mo) on the first factor P 2 {n,mo), we have 

dimT(/,re,mo) < dimQ{l,n,mo) — P 

< dimP 2 {n,mo) + nl — [mo/2] — P 

< dimP 2 {n,mo) — [mo/2] + l{n — 1) 

< dimP 2 (’^,”io) — [mo/2]-|-[re^/4] 

< dimP2(?^5^o) ~ 1 

where the last inequality follows from the definition of mo- This proves that T(/,re, mo) has 
codimension at least one in P{n, mo) for any 1 < / < re—1, therefore T(re, mo) has codimension 
at least one in P(re,mo) and this concludes the proof. □ 
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Proof of Provosition \A.3l To prove that the complement of SS{n,mo) in P 2 {n,mo) is con¬ 
tained in T(n,mo), we will prove that the complement of T(n, mo) in P 2 {n,mQ) is contained 
in SS{n,mo). So we fix Pq G P 2 {n,mo) \ T(n,mo) and 1 < / < n — 1. We denote by 0{n,l) 
the subset of Li{n,l) consisting of matrices whose columns are orthonormal vectors for the 
Euclidean scalar product. Recalling that the Grassmannian G{l,n) is the quotient of Li{n,l) 
by GL{1), it is also the quotient of 0(n,l) by the group 0{l) of orthogonal matrices of 
Therefore given any Aq G G{l,n), there exist an open neighborhood B\g of Aq in G{l,n) and 
a continuous map T : Raq 0{n,l) such that, if vr : 0{n,l) — ^ G{l,n) denotes the canonical 
projection, then tt o T is the identity. Let us now consider the continuous map 


F : P 2 {n,mo) x Raq ^ P 2 {l,mo), F{P,A) = P^(a)- 

Since Pq does not belong to T(n,mo), by definition of the latter set it comes that F{Po,A) 
does not belong to S(Z,mo) and therefore, by Theorem lA.ll F{Pq,A) G SE{l,mQ) for any 
A G i?Ao- Hence, by definition of SE{l,mo), there exist a neighborhood Ui of F{Pq,A) in 
P 2 {l,m) and positive constants C/,<5^ such that for any Q G Ui, the inequality 

max niin ||VQ(y)|| > GlC°~^ 

holds true for all 0 < ^ < 5^. Now by continuity of F, we can find a neighborhood V) of Pq 
in P 2 {n,mo) and an open neighborhood B'j^^ C Rao of Aq in G{l,n) such that F{Vi x Raq) 
contained in Z7;. So for any P G Vi and any A G Raq’ have 

inax niin ||VF(P, A)(y)|| > 

0<r?<S \\y\\=v 

for all 0 < .^ < d;. Now since the columns of the matrix T(A) form an orthonormal basis of 
A, setting x = T(A)y, x G A, ||x|| = ||y|| and hence 


min ||VF(P,A)(y)|l = min ||nAVP(x)ll = min ||VPA(a:)|| 

\\y\\=V l|a;||=T?, a:6A \\x\\=y,x&K 

where Ha is the orthogonal projection onto A, and Pa is the restriction of P to A. Therefore, 
for any P G V) and any A G Paq’ have 

max min ||VPa(x)|| > G[C°~^ 

0<ri<^ ||3;||=?7, xGA 


for all 0 < .^ < d;. To conclude, since the Grassmannian G{l,n) is compact, it can be covered 
by a finite number of neighborhoods of the form Pa^, Aq G G{1, n), and hence one can certainly 
find positive constants Gi,5i such that for any P G Vi and any A G G{l,n), the inequality 


max 

o<v<i 


min \\VPa{x)\\>GiC°~^ 

x\\=ri, IgA 


holds true for all 0 < ^ < 5 /. This means that Pq G SS{n, mo), and this finishes the proof. □ 


B Birkhoff normal forms with estimates 

The goal of this section is to give the proof of Proposition 12.51 using the work of Delshams 
and Gutifeez l |DG96j l. 
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Given I G N and P a homogeneous polynomial in ^ of degree I, if P{^) = Yl\u\=i ■: 
define the norm 

Ill’ll 

\v\=l 

By our analyticity assumption on the Hamiltonian H in (11.11) , we have the following expansion 
at the origin 

n 

= ^ HiiO = i ^ ^ Hi{^) 

l>2 j=l l>3 

and there exist positive constants c and d, which depends only on n, R and ||H^||r such that 
for any integer I >2, 

\\Hi\\ < c''-^d. (B.2) 

Using Cauchy formula one easily proves that 

\\Hi\\ < {2R)-^{e{2n + l)y\\H\\R 

and therefore one can choose 

c:= {2R)-^e{2n + l), d := {2R)-^{e{2n + l)f\\H\\R. (B.3) 

Given any function / that can be written as / = 'YhkPk^ with each P^ homogeneous of 
degree k in one easily check that 

sup 1/(01 (B.4) 

and, if g = fceven Qk-, with each Qk homogeneous of degree k/2 in /(Oi then 

sup \g{I)\ < \\Qk\\{p'^/2)’^/^. (B.5) 

k, k even 

Moreover, the above estimates hold true if / is replaced by a tensor-valued function. Recall 
the definition of 'I'q given in (II.5p . Recall that we also defined for any integer j > 3, ipa = 
0^=3'ha (i) and for convenience, we set := 1. We can finally state the main technical 
proposition of [DG96| . 

Proposition B.l (Delshams-Gutierrez). Let H be as in (jl.ip with a as in (|1.5I) and consider 
an integer K > A. If we define 

PK ■■= {MSncdKm{K))-^, 

then there exists a real-analytic symplectic transformation = Id + 0(0) defined on Bpj^ 
such that H o is in Birkhoff normal form up to order K, that is 

Ho^^{^) = a-m+ Y hk{m)+ Y 

k even, 4<k<K k>K-\-l 

where hk is a homogeneous polynomial of degree k/2 in I{fi), fk a homogeneous polynomial 
of degree k in f, with the following estimates: 

\\hk\\ < 6“^(6cd)^“^(/c — ^ even, A< k < K] 

WfkW < 20d^{20cd)^-^{K - 3)1{K - A: > A' + 1. 
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This is exactly the statement of Proposition 1 in [DG96] , to which we refer for the proof. 
We will now arrange these estimates in a way that will be more convenient for us. 

Proposition B.2. Let H be as in (HI) with a as in (irsi). Given an integer p > 2 and 
K > 2p, we have the following estimates on the homogeneous polynomials of Provosition 

\\hk\\ < 2p<k<K; (B.6) 

where 

(3{p) := 6-\6cd)^P-\2p - (B.7) 

and, given an integer 0 < q < K — 4, we have 

IIMI k>K + l. (B.8) 

where 

P{q) ■.= c-^d{2f)cdY{q + 2)\fjl+‘^. 

The proof of Proposition IB.21 is straightforward from Proposition IB. 11 Now from these 
estimates on the homogeneous parts of h™' and f^, we will deduce the estimates of Proposition 

ESI 

Proof of Proposition \2.5\. Recall that 

*”(/({)) = E 

k even, A<k<K 


where each hk is homogeneous of degree k/2. For p > 2 and k > 2p, V^hk is a tensor-valued 
homogeneous polynomial of degree {k — 2p)/2 and one can easily check (see |DG96j . estimates 
(24)), that 

iiv^h,ii < ikmhkW. 

Using this inequality, inequality (lB.5p and the estimate (IB.Op we get 


l|V^h-||, < 
< 

< 

< 

< 


\\vPhk\\{p^/2)^^-^py^ 

k even, 2p<k<K 

{k/2nhk\\{pV2)(>^-^py^ 

k even, 2p<k<K 

k even, 2p<k<K 

m Y {kl2Y{ll2)^^-^P)l\plpK)'^-^^ 

k even, 2p<k<K 

b{p) 


since the sum can be bounded by the corresponding series which is convergent. The same 
bound applies to ||V-^h”*||p for any j such that 2 < j < p, hence 

||V^/i”^||p_p = max WV^h^Wp < h{p). 
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Concerning 


f^{c)= E 

k>K+l 

since V/^ is a vector-valued homogeneous polynomial of degree A: — 1, we have 


l|V/fe||<A:||A|| 

and so, using this inequality together with inequality ()B.4I1 and the estimate (IB.Sp we obtain 


W^f^Wp 


< E I^ E 

k>K k>K 

< m E < mp^ E ipiPK?~‘'~^ 

k>K k>K 

< b{q)p<ie^. 


This concludes the proof. 


□ 


C Proof of the normal form statement 

C.l Technical estimates 

We first derive technical estimates for real-analytic vector fields defined on certain domains 
in C^"'. These estimates are stated and proved for Hamiltonian vector fields, even though the 
Hamiltonian character plays absolutely no role here. 

For 0 < j < n— 1, recall that ujj G M"'\{0} are Tj-periodic vectors, and that oj-i = 0 G M”. 
We write lwj{z) = uj ■ I(z), for 1 < j < n — 1, and we define the complex vector space 

Aj = {u G C” I u • uj-i = V ■ ujQ = ■■■= V ■ ojj-i = 0}. 

Then we consider three sequences of positive real numbers rj, and Sj, a sequence of points 
Zj G Brj and we let Xj = I{zj) + Aj be the complex affine subspace associated to Aj passing 
through I{zj). The complex domains we consider are given by 

= {^ e I Hz) G Aj, \I{z) - I{zj)\ < Sj, ||2:|| < Tj+^j}. 

We fix 0 < fjj < Sj and 0 < pj < ^j , and a real-analytic Hamiltonian vector field defined 
on Vsj,^Hzj)- Throughout this section, we will make the following two assumptions: 

CTj < Hj +Cj)pj, {L-i,Xj} = {L-oXj} = ■■■ = {Lj_i,Xj} = 0. (C.l) 

Lemma C.l. Assume that (jC.ip is satisfied. Then : Vsj-aj,^j-pjizj) '^sj,^j{zj) is a 
well-defined symplectic real-analytic embedding for all |t| < Tj = Hj + f,j)~^o'j\\X^.\\f^^^^, with 
the estimate \ \ - Id 11 ^< 11111 X^. \ \ . 

Proof. Let z G Vsj-aj,^j-pjizj) and z{t) = Xl^.{z) for small |t|, and let |s| < \t\. Since 
{Li,Xj}izis)) = ivi ■ {{h,Xj}izis)),... ,{In,Xj}izis))) := w/ • {/,Xj}(2(s)) 
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for — 1 < / < j — 1, the second part of (IC.ll) implies that 

■ {I,Xj}(.z{s)) = LOO ■ {I,Xj}{z{s)) = ■■■ = LOj-1 ■ {I,Xj}{z{s)) = 0, 
so {/, Xj}(-2('S)) G Aj which implies that 

d 


ds 


I{z{s)) = {I,Xj}{z(s)) G Aj 


and therefore 

I{z{t)) = I{z) + f ■^I{z{s))ds = I{z) - I{zj) + I{zj) + f ■^I{z{s))ds e Xj. 

Jo ds Jo ds 

Then, using the first part of (IC.lh . for 

|t| < (rj + = tj, 

we have 

\\z{t) - 2|| < \t\\\X^^\\rj,^j < Pj 
and, using Cauchy-Schwarz inequality, 

\I{z{t))-I{z)\ < 2-^\\z{t) + z\\\\z{t) - z\\ < {rj + ^j)\\z{t) - z\\ < irj+^j)\t\\\X^.\\s.,^^ <aj. 

This proves that X^. : Vsj-crj,^j-pjizj) —)• Vgj^^^izj) is a well-defined symplectic real-analytic 
embedding for \t\ < tj, with the estimate \\X^. - □ 

Lemma C.2. Assume that (jC.ip is satisfied, and let Xf be a real-analytic Hamiltonian vector 
field defined on Vsj^^jizj). Then, for |t| < rj/3 = (3(rj + ?j))~^<7j||-^xillrV’ have 


and therefore 




t 

Xj- 


Proof. Let \t\ < Tj/3 = (3(rj , We have the following expression 

(X^rx, = (DX-> o x‘,).(X^ o x‘,) = {dx-; o xl - Id) .{X, O X‘,) + X/ o X 
Lemma EJ] implies that X*^ : (aj) -t Vs^_2„^/3,s,-2p^/3(ap hence 

° K, - < l|OV7 - Id||,,_2a,/3,£,-2p3/3. 

We claim that 

W^^x^ -Id||s._2a^/3,^^-2p,/3 < 3(rj +f,j)(yj ^ll^x/ 
while obviously, using Lemma lC.il 

ll^x/ = W^-Xj -M|U,.-a^/3^-p,/3 < \t\\\X-x.\\sj,^j = 
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Assuming this claim, using the expression for {X^^)*Xf and putting all the estimates together, 
we arrive at 


^ (l + 3(r, +C,>. 


."I 


IX 


Xj 1i5j 


I ° JI -pj ’ 


therefore 




and also 

^ 2||X/l|s_2^./3,^._2p^Y3 

since |t| < rj/3 = (3(rj +^j))~^o'j\X^.\J^^^^. It remains to prove the claim. Let F = X~^ — Id, 
-2(Tj/3,^j-2pj/3(2j) aiid V G C^"’ a unit vector. The map 


eGC^F,,,(0=T(z + eu)GC2- 


is holomorphic for |^| < (3(rj +^ 7 )) Vj < pj/3, with + G ._^^Y 3 ^g_p./ 3 (z 7 ). The usual 
Cauchy’s estimate implies that 


||L)F( 2 ;)|| = sup \\DF{z).v 
lhll=i 


sup ||T^,J0)|| < 3(rjsup 
lhll=i l?l<( 3 hj+L)) 


hence 

|| DF (^) 1 | < 3 ( r , + e 7 ) u 7 '|| Tll ,^._,^/ 3,^^._,^./3 

and the claim follows since z G Vs^_2(7j/3,^j-2pj/3('2^i) was arbitrary. □ 

Lemma C.3. Assume that (jC.ip is satisfied, and letXf be a real-analytic Hamiltonian vector 
field defined on Vsj,(^jizj). Then 

11 > ^Xj ] 11 Sj -2(7j/3,^j-2pjl3 — 9 (’’i + ) Uj II ^Xj 11 Sj 11 11 Sj ■ 

Proof. We have the expression 

= fAx^rx, 

t=0 

so for z G Vsj-(Tj,^j-pjizj), let us define the holomorphic map 

t G C ^ F,{t) = (X* .)*X/(z) G 

for |t| < Tj/3 = {3{rj + ?j))~^Uj||Xj^^. . By Cauchy’s estimate 

||[Xy,X;,J(z)|| = ||Fi(0)|| < 3r-i||F,(t)|| < 
and by Lemma IC.21 

3r-i||(X‘^.)*X;||,^._.^,^^_,^, < 6r-i||X;||,^._2.^/3,^,-2p,/3 < 6(r,-+C7)u7i||A^J|,^,,^J|X;||,^,^^. 
Since z G Vsj-aj,^j-pj{zj) was arbitrary, this proves that 

\\[Xf,Xx,j]\\sj-(7j,^j-Pj < 6(rj +Cj)Uj ^||X;^Jls^^^J|X/||s^.^^^. 

and the lemma follows by simply replacing aj and pj by respectively 2ajj3 and 2pjl3. □ 
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Lemma C.4. Assume that (jC.ll) is satisfied, and let he a real-analytic Hamiltonian vector 
field defined on Vrj^^j{zj), which is integrahle, that is k is a function of I{z) alone. Then 

||[^A:)^Xj]llsj-2(Tj/3,^j-2pj/3 — ^Pj I I^Xj I Isj 11 I Isj ' 

Proof. Here we write 

t=o 

and we observe that since X/. is integrable, if we let z(t) = Xj:(z), then I(z{t)) = I{z). 
This implies that Xj, : 'Ds.^^._p.{zj) —)• Vg.^^.^Zj) is a well-defined symplectic real-analytic 
embedding for all \t\ < r' = pj\\X]f\\~^c.. The conclusion follows easily by repeating all the 
previous arguments with rj instead of tj. □ 


C.2 Proof of the Proposition 13.11 

Proposition 13.II will be proved by iterating m times an averaging procedure, which is classical 
in the case j = 0 , but more involved in the general case. 

Proof of Provosition fOl Let us fix 0 < j < n — 1, and set £j := 2^e. The integer m > 1 
being given, for 0 < i < m we define 

e* := 7 * := (1 - 2 “*) 2 ej, s't := 3sj - isj/m, 

Then we claim that for each 0 < i < m, there exists a real-analytic symplectic embedding 
'^ 3 sj, 3 ^j{z:j) such that 

{Hj- ffio^) = {h + gfio^) = h + gi+fj 


with 


,g^} = {Lo ,9j} = --- = {L ,, 4 } = 0 , {L _,, f^} = {L ,, /*} 

and with the estimates 


{L. 


-1 5 


i^} = o 




s’-. 

3 ’^3 


< T 


\X 


■ ■ < P 
s),£] - ^3’ 


I4>* 
I J 


-Id||, 


.1 Cl 

'3’P 


<T3l3- 


Let us prove the claim by induction on 0 < z < m. 

For z = 0, letting be the identity, s'? := 0 and /? ■= Qj, there is nothing to prove. Then 
assume that the statement holds true for some 0 < z < m — 1, and let Ffj = {Hj — fj) o = 
h + Oj + fj ■ We define the functions 


■= T, 


-1 


noXfdt, ^:=T: 


'0 Jo 

whose associated Hamiltonian vector fields are given by 


^ £'t{rj-[fj]3)°K3dt 




fT 

Jo 


*ndt, x.=T-^ 


t{xiy{n-[f]]j)dt 
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with, using our inductive assumption, the following obvious estimates 




It is clear that 


{[/;].,m=o. 


(C.2) 

(C.3) 


For — 1 < / < j — 1, observe that L, } = 0, so o X^j. = and hence 


[/*],} = 


Jo 


fioXlMt 


= 


-1 


10 


{L,oXl^,noXl\dt 


= T-^ l^\L,J}}oXiM 


where the last equality follows from the symplectic character of X^^.. Using our inductive 
assumption, this implies 


{L_1, [/]],} = {Lo, im = ■■■ = {4-1, [/llil = 0, 

and by a completely similar argument, we also get 

{L_i,x4 = {4,Xi} = ••• = = 0. 

Now set 


(C.4) 


(C.5) 


aj := Sj -sy^ = Sj/m, pj := Q - = ij/m. 


Since , using the first inequality of (13.3p we have 

(Tj = Sjjm < 4 + 2^j)^jjm < {rj + Cj)Cj/m = (rj + Cj)Pj 
and therefore, using also (|C.5p . we can apply Lemma fC.II 




is a well-defined symplectic real-analytic embedding for all 

|i| < Tj = 4 = {m{rj + Cj))~^Sj\\X^^\\-^ i, 

3 


with the estimate llXh — Idll i+i ,i+i < |t|||X i || i «. Moreover, as Ci < 3^,-, using the second 
_ Xj U ’U U j 

estimate of (IC.2p and the second inequality of (13.311 . we have 


Tj > (m(rj-|-3^j)Tje*) = 2^{m{rj+‘i^j)Tj2^e) ^Sj > {m{rj+2>^j)Tj2^£) > 216 (C.6) 


so Tj > 1 and hence Xh : V^i+i y+i{zj) —)• Vgi^^iXzj) is well-defined, with 


(C.7) 
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It is easy to check, using an integration by parts, that {xj, luij} = /j — [/jljj and this equality, 
together with Taylor’s formula with integral remainder gives 

{h + gj + fj) o X^i = h + Qj + [fj]j + fj 


with 




We set o Xh, 5*+^ = g) + [f-]j and = /) so that 


{Hj -fj)o d-f 1 = m o = {h + g] + n) o xh^ =h + 5!+^ + /: 


,2+1 I ^2 + 1 




First 


observe that +• : V^^i+i ^i+i{zj) —)• 'D^sj,2,ij{zj) is a real-analytic symplectic embedding, 

and using ()C. 7 p together with our inductive assumption, we have the estimate 


1$ 


i+l 


- Id|L.+i,g.+i < 11+ -+ ||Xb -Id||^.+I_g ,+1 < TjiYj + e]) = r,7;+i 


Then 

= {+++'} = • • • = {+++'} = 0 

follows from the definition of <7+^, the inductive assumption, (IC.3D and (IC.4p . Moreover, 
using the first estimate of (IC.2h and our inductive assumption, 

For — 1 < / < j — 1, we already know that 

{++} = {+, f^} = {L,, [f^W = {+ ,x)} = o 

which implies = 0 whereas {luii,h — 1^^} = 0 is obvious. These equalities, together 

with Jacobi identity, imply that 


{+ ,{{h - +.) -I- 5 * -I- fjj-, Xj}} — 0 


and therefore 


{+,/+'}= / {L„{{h-L,) + g} + flt,x]}}oXl,dt. 


It follows that 

{c_., /]+'} = {;„.,/■+') = ■ ■ ■ = {c,.,, /■+') = 0 . 

To complete the proof of the claim, it remains to estimate 


X ri+l — 
■’3 


[\x^,r[Xh-i^ + Xg. + Xf. X^.]dt. 


First, 




r+l|L*+l £»+l < sup ||(X*i)*[X/i_;,^ -h X i + X,i X i]|| i+l ,i+i 

Jj 7 Xj 3 yj Xj 7 Aj 
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and since r > 3 by ()C. 6 |) . we can apply Lemma [0.21 to get 


+"^95 + ^fi^^xi]\\s)-2aj/3,G-2pj/3- 


Using Lemma IC.31 we have 


and since 
we get 


11 fit ’ ] 1 1^5 /U] -2pi /3 ^ 9 (^i + €' ) O-j II ^fi^ 1 1 j 11 11 s} 4 


ll[^/L’^x5]ll^1-2-7,/3,q-2p,/3 < {9{rj +3Cj)mTjejSj 


-^j j 31 ^3 

since < 3^j and e*- < Sj. Similarly, since 7 *- < 2ej, 


(C. 8 ) 


(C.9) 


ll[^q.>^x']ll^'- 2 cT,/ 3 ,e- 2 ft /3 < {9{rj +3ij)mTj-i]Sj )£) < {lS{rj + 3ij)mTj£jSj )e*. (C.IO) 


'■i ' '^’‘>3 ■ 

Concerning the last bracket, let us first prove that 


[X^.i X^.] = [Xh, -I X ] 

3 '^3 3 

where we recall that Xj = I{zj) + Kj and h^, is the restriction of h to Xj. To do this, it is 
sufficient to prove {h — lu]j,x''j] = {h\. — LjjXj}) which is equivalent to {/i,y*} = 

For any z £ V^i we have 

{h,4}(z) = vh(/(z))-{/,x;.}(z). 

But for any — 1 < f < j — 1, we know that 


{^ioi,x]]{z) = oJi ■ {I,x]]{z) = 9 

which means that {I,x^j}{^) ^ ^j- Therefore, recalling that ITj denotes the orthogonal 
projection onto Aj, it comes that 

{Kxi]{z) = %{Vh{I{z))) • {I,x)]{z) = Vh~^^{I{z)) • {I,x]]{z) = {h-^.xi]{z) 

and therefore {h,x)} = {^x-iX]}- Then, for any z £ we can estimate 


\\Xh-^^{I{z))-uj,\\ = 

< 

< 

< 

< 


\\tij{Xh{I{z))) - ojjW 

\\ii,{xh{i{z))) - ii,{Xh{i{zm\ + l|n,(vM/(^i))) 

\\Xh{I{z))-Xh{I{zj))\\ + \\UjiVh{I{z,)))-u;j\\ 
F'WHz) - ^( 2 : 9)11 + Sj < F^/n\I{z) - I{zj)\ + Sj 
Fy/ns^j + Sj < {3Fy/n + l)sj 


UJi 


since s*- < 3sj, and where we used the fact that 


sup \\X^h{I{z))\\<F 

(Zj) 
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From this, we deduce that 


-^‘^j)iHz))\\\z\\\ < {^F^/^+l)sj{rj+ Cj) 

' " ^ ^ zev^i i {zj) 

j j 

< {3F^/n + l)sj{rj + 3^j) 

and using Lemma IC.4[ we get 

11 , X^]\\si-2aj/3,ej-2pj/3 = 11 , X^^\\si-2aj/3,ej-2pj/5 

< 9pj 11 X^i 11 

hence 

\\[^h-iu,.,Xx*^^^s^j-2crj/3,Cj-2pj/3 ^ ^Pj ^{3F^/n + l)sj{rj + 3ij)Tje^j 

= {9{3Fy/n + + 3(j)mTjSj)e'j. (C.ll) 

Putting the estimates (|C.8p . (IC.9p . PC.IOI) and pC.lip together, and recalling that £j = 2%, 
we arrive at 


(2-^54(rj + 3^j)mTj£Sj ^ + 18{3F^/n + ^{rj + 3^j)mTjSj)£'j. 

Using the second and third inequality of (13.3p . we obtain 

I |^f*+i I|g»+i g»+i ^ . 

Jj j ’’ j '' 

This finishes the proof of the claim. 

Now let us define g~^ = gF and = f^ + fjo^j. Since s™ = 2sj and = 2^j, 

is a real-analytic symplectic embedding 

■ V2sj,2^j{Zj) V3sj,3^j{Zj) 

such that Hj o = h + g'^ + f^. We already know that {/oj-dS'^} = = ■■■ = 

1 = 0) the estimates 

\\Xg+\\2s„2i, <lT< 2ej = 2^-+ie, ||4>o - Id||2.„2€, < ^^'7™ < 2^’+'T;-e. 

To conclude the proof of the proposition, it remains to estimate Xf+. First recall that 

/ 7 


\\XfY' 


.<£? = 2-™e, = 2^2- 




e. 


(C.12) 


Then, fj o = fj o 4>"* = fj o X\ o • • • X^m and so Xf.o,i>^ = (X^m)* ■ ■ ■ {X\)*Xf ., where 



MIL 


i + 1 ci + 1 

'j 


<11^. 


I ■ ■ < T-p® — 


2“®T f • 

Z, J-jtj, 
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for each \ < i < m. Applying Lemma IC.21 inductively yields 


IX 


fjo^j I \ 2sj ,2^j 


< 

< 

< 


m—1 


s n (i + 3 ('-i + Q)y'll^x‘ll.;.q) Ui,h.,a(, 

2=0 

m—1 

(^1 + 3(rj + I 

i=o ^ 

m—1 

n (1 + 2 *3(rj + 2>ij)s- ^'mTjej)\\Xf^\\2s^,2ij 




i=0 


^ m—1 


exp I 2 *3(rj + ^Cj)Sj j ||X 


fj \ \ 2sj,2^j 


i=0 


exp(6(rj + 3^j)sj mTjej)\\Xfj \ | 2 s^., 2 C, 


The second condition of (13.31) implies in particular that exp(6(rj + 3^j)sj ^rriTjEj) < 2 and 
therefore 




\2s,,2i, < 2||XyJ|2.,,2^, < 2(j2^-i)2-™e = j2^2-^e. 


(C.13) 


From (jC.12h and ()C.13p we get 


||X^+||2.,,2?, < 2^2-”*e +j2^2-™e = (j + l)2^2-™e, 
and this ends the proof. 


□ 
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